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' In quantum field theory, there exist states in which the expectation value of the energy density for 

a quantized field is negative. These negative energy densities lead to several problems such as the 
Q"^ ! failure of the classical energy conditions, the production of closed timelike curves and faster than 

I light travel, violations of the second law of thermodynamics, and the possible production of naked 

O^' singularities. 

^ , Although quantum field theory introduces negative energies, it also provides constraints in the 

form of quantum inequalities (QI's). These uncertainty principle- type relations limit the magnitude 
and duration of any negative energy. We derive a general form of the QI on the energy density for 
^ . both the quantized scalar and electromagnetic fields in static curved spacetimes. In the case of the 

H ' scalar field, the QI can be written as the Euclidean wave operator acting on the Euclidean Green's 

function. Additionally, a small distance expansion on the Green's function is used to derive the QI 
in the short sampling time limit. It is found that the QI in this limit reduces to the flat space form 
with subdominant correction terms which depend on the spacetime geometry. 

Several example spacetimes are studied in which exact forms of the QI's can be found. These 
include the three- and four-dimensional static Robertson- Walker spacetimes, fiat space with per- 
fectly refiecting mirrors, Rindler and static de Sitter space, and the spacetime outside a black hole. 
In all of the above cases, we find that the quantum inequalities give a lower limit on how much 
negative energy may be observed relative to the vacuum energy density of the spacetime. For the 
particular case of the black hole, it is found that the quantum inequality on the energy density is 
measured relative to the Boulware vacuum. 

Finally, the application of the quantum inequalities to the Alcubierre warp drive spacetime 
leads to strict constraints on the thickness of the negative energy region needed to maintain the 
warp drive. Under these constraints, we discover that the total negative energy required exceeds 
the total mass of the visible universe by a hundred billion times, making warp drive an impractical 
form of transportation. 
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Chapter 1 

Energy Conditions in General 
Relativity 

1.1 The Classical Energy Conditions 

The most powerful and intriguing tool in cosmology is Einstein's equation, 

G^v = R^j,u — -Rg/iu = SirGT^iy. (1.1) 

Here we have chosen the spacetime metric to have a (—,+,+,+) signature and the definitions of 
the Riemann and Ricci tensors in the convention of Misner, Thorne and Wheeler Einstein's 
equation relates the spacetime geometry on the left with the matter source terms on the right. 
The equation is interesting in that it admits a large number of solutions. If there are no physical 
constraints on the metric or the stress-tensor, then Einstein's equation is vacuous since any metric 
would be a "solution" corresponding to some distribution of the stress-energy. We could then 
have "designer spacetimes" which could exhibit any behavior that one likes. Therefore, if we can 
place any limitations upon the equation, or on terms in the equation, we could limit the search for 
physically realizable solutions. 

The most trivial condition comes from the Bianchi identities, that the covariant divergence of 
the Einstein tensor must vanish: 

Vf,G^"' = 0. (1.2) 
From Einstein's equation, this leads to the condition on the stress-energy tensor of 

This is the covariant form of the conservation of energy in general relativity. 

The next condition which is generally considered reasonable for all classical matter is that the 
energy density seen by an observer should be non-negative. Consider an observer who moves along 
the geodesic 7^(r). Then his four- velocity, 

^M(^) ^ ^jf^ir), (1.4) 

is everywhere timelike and tangent to the curve 7^(r). In the observer's frame, the statement that 
the energy density should be non-negative is 

Tf,^u''{T)u''{T) >0, for all timelike ?i''(T). (1.5) 
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This is typically called the Weak Energy Condition (WEC). This condition was an important 
assumption in the proof of the singularity theorem developed by Penrose In the early part 
of this century, it was known that a spherically symmetric object that underwent gravitational 
collapse would form a black hole, with a singularity at the center. However, relativists were unsure 
if the same would be true for a more generic distribution of matter. Penrose's original singularity 
theorem showed that the end product of gravitational collapse would always be a singularity without 
assuming any special symmetries of the matter distribution. The WEC is used in the proof of 
the singularity theorem in the Raychaudhuri equation to ensure focusing of a congruence of null 
geodesies. An additional assumption sometimes used in the proof of the singularity theorems is the 
Strong Energy Condition (SEC), 

(Tf,u - ^^5/^^^) u^'{t)u''{t) > 0, for all timelike u^ir). (1.6) 

Here T = Tfj_^ is the contraction of the stress-tensor. If the SEC is satisfied, then it ensures that 
gravity is always an attractive force. The SEC is "stronger" than the WEC only in the sense that 
it is a more restrictive condition, and therefore it is much easier to violate the SEC than the WEC. 

Immediately following the successes of Penrose's singularity theorem for black holes. Hawking 
1^, ^ developed a similar theorem for the open Friedman-Robertson- Walker universe which showed 
that there had to exist an essential, initial singularity from which the universe was born. Nearly 
thirty years later, singularity theorems are still being proven for various spacetimes. For example, 
Borde and Vilenkin have proven the existence of initial singularities for inflationary cosmologies 
and for certain classes of closed universes |5|, ^ . 

For completeness, we also mention two other classical energy conditions. The first is the Null 
Energy Condition (NEC), 

T^^Rf^K" > 0, for all K" , (1.7) 

where is the tangent to a null curve. The NEC follows by continuity from both the WEC ( |l.5D 
and the SEC ( |1.6| ) in the limit when becomes a null vector. The second condition is a little 
different from the previous ones. Let us assume, as we did above, that is a future-directed 
timelike vector. It follows that the product, —T'^^u'^, should also be a future-directed timelike or 
null vector. This is known as the Dominant Energy Condition (DEC) and can be interpreted to 
mean that the speed of energy flow of matter should always be less than the speed of light. 

We see that the energy conditions at the classical level are an extremely useful tool that have 
yielded some remarkable results in the form of the singularity theorems and positive mass theorems. 
However, many of the energy conditions are doomed to failure when quantum field theoretic matter 
is introduced as the source of gravity. A first approximation to a "quantum theory of gravity" is 
achieved via the semiclassical equation 

G^, = 8ttG{T^,). (1.8) 

We are still treating the spacetime as a "smooth" classical background, but we have replaced the 
classical stress-tensor with its renormalized quantum expectation value. It is presumed that the 
field equation should be valid in the test field limit where the mass and/or the energy of the 
quantized field is not so large as to cause significant back-reaction on the spacetime. However, 
the replacement of the classical stress-tensor with its quantum expectation value is not without 
difficulty. At the same time that Penrose and Hawking were introducing the singularity theorems, 
Epstein, Glaser and Jaffe were demonstrating that the local energy density in a quantized field 
theory was not always positive definite. We will see in Chapter § that it is possible to construct 
specific quantum states in which the energy density is negative along some part of an observer's 
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worldline. The simplest example, discussed in Section 2.1, is a state in which the particle content 
of the theory is a superposition of the vacuum plus two particle state. Another example is the 
squeezed states of light which have been produced and observed experimentally although it is 
unclear at this time whether it is possible to measure the negative energy densities that exist for 
these two cases. However, there is the interesting case of the Casimir vacuum energy density for 
the quantized electromagnetic field between two perfectly conducting parallel plates [^. Here the 
vacuum energy between the plates is found to be a constant and everywhere negative. Indirect 
effects of the negative vacuum energy, in the form of the Casimir force which causes the plates to 
be attracted to one another, has been measured experimentally 11]. It has even been suggested 



that the Casimir force may play a role in the measurement limitations in electron force microscopy 



|2|. 



Related to the Casimir vacuum energy is the vacuum energy of a quantized field in curved 
spacetime, which we will discuss further in Section An example is the static closed universe, 
where the vacuum energy for the conformally coupled field is positive p^ , while that of the 
minimally coupled field is negative |]l5|. Instances of non-zero vacuum energies in the semiclassical 
theory of gravity seem to be extremely widespread. The most well known, and probably most 
carefully studied, is the vacuum energy around a black hole. Hawking has shown that this causes 
the black hole to evaporate, with a flux of positive energy particles being radiated outward. In this 
case both the WEC and the NEC are violated ||, 0, |l|, ||, || . 



In terms of a consistent mathematical theory of semiclassical gravity, the existence of negative 
energy densities is not necessarily a problem. On the contrary, it could be viewed as beneficial 
because it would allow spacetimes in which time travel would be possible, or where mankind 
could construct spacecraft which could warp spacetime around them in order to travel around the 
universe faster than the speed of light. Even wormholes would be allowed as self-consistent solutions 
to Einstein's equation. 

Objections to negative energies are raised for several reasons such as violations of the second law 
of thermodynamics, the possibility of creating naked singularities by violating cosmic censorship, 
and the loss of the positive mass theorems. On the more philosophical side, if time travel, space or 
time warps, wormholes and any other exotic solution that physicists can dream up were allowed, 
then we would have to deal with a new set of complexities, such as problems with causality and 
the notion of free will. Visser's description of four possible schools of thought on how to build 
an entirely self-consistent theory of physics with regard to time travel is quoted here: 

1. The radical rewrite conjecture: One might make one's peace with the notion 
of time travel and proceed to rewrite all of physics (and logic) from the ground up. 
This is a very painful procedure, not to be undertaken lightly. 

2. Novikov's consistency conjecture: This is a slightly more modest way of mak- 
ing one's peace with the notion of time travel. One simply asserts that the universe 
is consistent, so that whatever temporal transpositions and trips one undertakes, 
events must conspire in such a way that the overall result is consistent. (A more 
aggressive version of the consistency conjecture attempts to derive this "principle 
of self-consistency" from some appropriate micro-physical assumptions.) 

3. Hawking's chronology protection conjecture: Hawking has conjectured that 
the cosmos works in such a way that time travel is completely and utterly forbidden. 
Loosely speaking, "thou shalt not travel in time," or "suffer not a time machine 
exist." The chronology protection conjecture permits spacewarps/wormholes but 
forbids timewarps/time machines. 
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4. The boring physics conjecture: This conjecture states, roughly: "A pox upon 
aU nonstandard speculative physics. There are no wormholes and/or spacewarps. 
There are no time machines/timewarps. Stop speculating. Get back to something 
we know something about." 

Besides the philosophical reasons, there is one particularly important mathematical result of 
the classical theory of gravity that would fail if negative energies did exist and were widespread. 
These are the singularity theorems which establish the existence of singularities in general relativity, 
particularly for gravitational collapse and at the initial "big bang." 



1.2 The Averaged Energy Conditions 

In the 1970's and 1980's, several approaches were proposed to study the extent to which quantum 
fields may violate the weak energy condition, primarily with the hope of restoring some, if not all, of 
the classical general relativity results. In 1978, Tipler suggested averaging the pointwise conditions 
over an observer's geodesic Subsequent work by various authors |2^, ||, |2|, |2|, pOl p| 



eventually led to what is presently called the Averaged Weak Energy Condition (AWEC), 

{Tf,^u''u'')dT > 0, (1. 



where is the tangent to a timelike geodesic and r is the observer's proper time. Here the energy 
density of the observer is not measured at a single point, but is averaged over the entire observer's 
worldline. Unlike the WEC, the AWEC allows for the existence of negative energy densities so long 
as there is compensating positive energy elsewhere along the observer's worldline. 
There is also an Averaged Null Energy Condition (ANEC), 



{T^^K''K'')d\ > 0, (1.10) 
where K^^ is the tangent to a null geodesic, and A is the affine parameter along that null curve. 



Roman [|2^, has demonstrated in a series of papers that the singularity theorem developed 
by Penrose will still hold if the weak energy condition is replaced by the averaged null energy 
condition. Therefore, within the limits of the semiclassical theory there exists the possibility of 
recovering some, if not all, of the classical singularity theorems. 

The ANEC is somewhat unique in that it follows directly from quantum field theory. Wald 



and Yurtsever |31| showed that the ANEC could be proven to hold in two dimensions for any 
Hadamard state along a complete, achronal null geodesic. They also showed that the same result 
held in four-dimensional Minkowski spacetime as well. The results in two-dimensional fiat spacetime 



were confirmed by way of an entirely different method by Ford and Roman [32|. The question of 
whether back-reaction would enforce the averaged null energy condition in semiclassical gravity has 
been addressed recently by Flanagan and Wald p^ ] . Their work seems to indicate that if violations 
of the ANEC do occur for self-consistent solutions of the semiclassical equations, then they cannot 
be macroscopic, but are limited to the Planck scale. On such a tiny scale, we do not know if the 
semiclassical theory is reliable. It is expected that at the Planck scale the effects of a full quantum 
theory of gravity would be important. 

At first the averaged energy conditions seem like a vast improvement. If the averaged energy 
conditions hold in general spacetimes, there is hope that we can restore the singularity theorems 
and gain a new class of energy conditions better suited for dealing with quantum fields. However, 
it is straightforward to show that both the AWEC and the ANEC are violated for any observer 
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when the quantum field is in the vacuum state in certain spacetimes. The simplest example is an 
observer who sits at rest between two uncharged, perfectly conducting plates. Casimir |^ showed 
for the quantized electromagnetic field that the vacuum energy between the plates is negative. 
Therefore, over any time interval which we average, the AWEC will fail. Similar effects can occur 
in curved spacetime for the vacuum state. Examples are the various vacuum states outside a black 



hole [O, IS, 19, 20, 21, 22|. It appears that a difference between an arbitrary state and the vacuum 



state often obeys a type of AWEC, which can be written as 

/oo 
[(VlT^.n'^nnV) - (0,|T^,^x^nnOe)] dr > 0. (1.11) 
-oo 

This equation, which has been called the Quantum Averaged Weak Energy Condition (QAWEC) is 
derived in Section \i.2\ It holds for a certain class of observers in any static spacetime |3^, |3^, ^ 
and is the result of the infinite sampling time limit of the quantum inequalities, which will be defined 
below. 

Another aspect of the averaged energy conditions is that there are a number of "designer" 
spacetimes where they are necessarily violated. Three examples of particular interest are: the 
traversable wormhole [p^, |3H|, the Alcubierre "warp drive" [p^, HH], and the Krasnikov tube |41, 42 



All three require negative energy densities to maintain the spacetime. Likewise, all three have 
similar problems in that slight modifications of these spacetimes would lead to closed timelike curves, 
and would allow backward time travel Thus, these spacetimes may be entirely unphysical in 
the first place. It has also been shown by various researchers that if warp drives and wormholes 
exist, then they must be limited to the Planck size ^2], 44 1; otherwise extreme differences in 



the scales which characterize them would exist. In Chapter ^ we will discuss the example of the 
macroscopic warp drive, where we find that the warp bubble's radius can be on the order of tens 
or hundreds of meters, but the bubble's walls are constrained to be near the Planck scale. Similar 
results have been found for wormholes and the Krasnikov tube. 

There is also a limitation in the averaged energy conditions. The AWEC and the ANEC say 
nothing about the distribution or magnitude of the negative energy that may exist in a spacetime. 
It is possible to have extremely large quantities of negative energy over a very large region, but 
the AWEC would still hold if the observer's worldline passes through a region of compensating 
positive energy. If we can have large negative energies for an arbitrary time, it would be possible to 
violate the second law of thermodynamics or causality over a macroscopic region Also, within 
the negative energy regions the stress-tensor has been shown to have quantum fluctuations on the 
order of the magnitude of the negative energy density p6[ | . This could mean that the metric of the 
spacetime in the negative energy region would also be fluctuating, and we may lose the notion of 
a "smooth" classical background. 



1.3 The Quantum InequaUties 

A new set of energy constraints were pioneered by Ford in the late 1970's [^], eventually leading to 
constraints on negative energy fluxes in 1991 ||4^. He derived, directly from quantum field theory, 
a strict constraint on the magnitude of the negative energy flux that an observer might measure. 
It was shown that if one does not average over the entire worldline of the observer as is the case 
with the AWEC, but weights the integral with a sampling function of characteristic width tQ, then 
the observer could at most see negative energy fluxes, F^, in four-dimensional Minkowski space 
bounded below by 

to r JTH^^ > 3_ 



F^ = - ^T-^dt > (1.12) 
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Roughly speaking, if a negative energy flux exists for a time to, then its magnitude must be less 
than about . We see that the specific sampling function chosen in this case, tQ/'K{t^ + tg), is 
a Lorentzian with a characteristic width to ai^d 8' time integral of unity. Ford and Roman also 
probed the possibility of creating naked singularities by injecting a flux of negative energy into a 
maximally charged Reissner-Nordstrom black hole ||l^. The idea was to use the negative energy 
flux to reduce the black hole's mass below the magnitude of the charge. They found that if the 
black hole's mass was reduced by an amount AM, then the naked singularity produced could only 
exist for a time AT < |AM|-^. 

Ford and Roman introduced the Quantum Inequalities (QI's) on the energy density in 1995 



|32|. These uncertainty principle- type relations constrain the magnitude and duration of negative 
energy densities in the same manner as negative energy fluxes. Originally proven in four-dimensional 
Minkowski space, the quantum inequality for free, quantized, massless scalar fields can be written 
in its covariant form as 

P = — — 1- 9- > - oo 2 4 ' "^0- (1-13) 



vrj-oo r2+r2 " 327r2r4' 

Here is the tangent to a geodesic observer's worldline and r is the observer's proper time. The 
expectation value () is taken with respect to an arbitrary state 1-0), and tq is the characteristic 
width of the Lorentzian sampling function. Such inequalities limit the magnitude of the negative 
energy violations and the time for which they are allowed to exist. In addition, such QI relations 
reduce to the usual AWEC type conditions in the infinite sampling time limit. 

Flat space quantum inequality- type relations of this form were first applied to curved spacetime 



by Ford and Roman [44| in a paper which discusses Morris-Thorne wormhole geometries. It was 
argued that if the sampling time was kept small, then the spacetime region over which an observer 
samples the energy density would be approximately flat and the flat space quantum inequality 
should hold. Under this restriction they showed that wormholes are either on the order of a few 
thousand Planck lengths in size or there is a great disparity between the length scales that char- 
acterize the wormhole. Inherent in this argument is what should characterize a "small" sampling 
time? In their case, they restricted the sampling time to be less than the minimal "characteristic" 
curvature radius of the geometry and/or the proper distance to any boundaries. When a small 
sampling time expansion of the quantum inequality is performed in a general curved spacetime, 
such an assumption is indeed born out mathematically. This short sampling time technique is 
particularly useful in nonstatic spacetimes where exact quantum inequalities have yet to be proven 
exactly. 



1.4 Curved Spacetime Quantum Inequalities 

Although the method of small sampling times is useful, it does not address the question of how the 
curvature will enter into the quantum inequalities for arbitrarily long sampling times. This is the 
main thrust behind the present work. In Chapter ^ we will derive, for an arbitrarily curved static 
spacetime, a formulation for finding exact quantum inequalities which will be valid for all sampling 
times. For the scalar field, we find that the quantum inequality can be written in one of two forms, 
either as a sum of mode functions defined with positive frequency on the curved background, 

Ap^*-^r ^i|^dt>-Ef^^ + >V,VA(x)|V^-\ (1.14) 

71 J— oo ~r "-O X \\ytt\ ^ J 
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or as the Euclidean wave operator acting on the Euchdean Green's function (two-point functions) 

A/5> -^□irGij(x,-to;x,+to). (1.15) 

Both forms are functionally equivalent, but the Euclidean Green's function method is often much 
easier to use if the Feynman or Wightman Green's function is already known in a given spacetime. 
If the Green's function is inserted into the expression above, we can find the general form of the 
quantum inequality. In three dimensions it is given by 

A/5 > --^^ (^)5(m, To,g^p) , (1.16) 

and in four dimensions. 



167rTQ 



A/5 > -iTTT^ ^'^5(m, TO, gap) ■ (1.17) 

Here, tq = I^mI^^^ to is the proper sampling time of a static observer, g^/s is the metric for the static 
spacetime, m is the mass of the scalar particle, and ^^^S and ^^^S are called the "scale" functions. 
They determine how the flat space quantum inequalities are modifled in curved spacetimes and/or 
for massive particles. 



In Section 3.2 we will show that in the inflnite sampling time limit, the quantum inequality will 



in general reduce to the condition 

hm — / ^ , 2 ^ Pvacuum- (1.18) 

Here we are taking the expectation value of the renormalized energy density on the left-hand side 
with respect to some arbitrary particle state and Pvacuum is the energy density of the vacuum 
defined by the timelike Killing vector. This is a modification of the classical AWEC inequality by 
the addition of the vacuum energy term. 



In Section |3.3| we will proceed with a small sampling time expansion of the Green's function to 
obtain an asymptotic form for the quantum inequalities. In this limit we find that the expansion 
of the quantum inequality reduces to the fiat space form, with subdominant terms dependent upon 
the curvature. 



In Section 3.4, we will develop the quantum inequality for the quantized electromagnetic field in 
a curved spacetime. The flat space quantum inequality was originally derived by Ford and Roman 
and takes an especially simple form 

The difference of a factor of two between the scalar field and electromagnetic field quantum inequal- 
ities occurs because the photons have two spin degrees of freedom, while the scalar particles have 
only one. As we will see, the electromagnetic field quantum inequality in static curved spacetimes 
can also be written as a mode sum, very similar to the mode sum form of the quantum inequality 
for the scalar field. At present, it is not known if there also exists a Green's function approach for 
the formulation of the quantum inequalities for the quantized electromagnetic field. 

In Chapter |^ we will apply the quantum inequality to several examples of curved spacetimes in 
two, three, and four dimensions. In two dimensions, the quantum inequality on the scalar field for 
static observers will be shown to have the general form 
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Again, tq is the proper sampling time of the observer. This is the inequahty in all two-dimensional 



spacetimes. In Section 4.1 we will see that this occurs because of the conformal equivalence between 



all two-dimensional spacetimes. 



In Section 4.2 we will find the exact quantum inequalities for the three-dimensional equivalent 
of the closed and flat static Robertson- Walker spacetimes. In this way, the explicit form for the 
scale function in Eq. ( |1.16 ) will be determined, and we will examine its behavior. This will be 



followed by the application of the quantum inequalities in four-dimensional static Robertson- Walker 
spacetimes. In all of the above cases, the process of renormalization will be important for the 
determination of the renormalized quantum inequality. In four-dimensional static spacetimes, such 
as those mentioned above, the renormalized quantum inequality is given by 

to {Ttt/ gtt) Ren. ^ 3 (4) , , , 

— / ,2 , ,2 - " on 2 4 '^i/^' ^0, ^a/?) + Pvacuum- (1-21) 

J — oo t + O^TT Tq 

Here the vacuum energy, Pvacuurm is important in determining what the absolute lower bound may 
be on any measurement of the renormalized energy density. 



In Section 4.4 we will consider the particular case of the quantum inequalities in flat spacetimes 
in which there exist perfectly reflecting mirrors. Both for a single mirror and for parallel plates, we 
will find that the quantum inequality is modified from the flat space form due to boundary effects. 
Not only is the spacetime curvature important in determining the form of the quantum inequality, 
but so is the presence of a boundary. We will continue along these lines by looking at spacetimes 



in which there are horizons. In Section we will look at quantum inequalities in Rindler and 
de Sitter spacetimes. 



In Section 4.6 quantum inequalities will be found for two- and four-dimensional black holes. 



For the black hole in two dimensions, we find the exact renormalized quantum inequality to be 

- "87r(l-2M/r)t2 + PBoulware- (1-22) 

The first term is the standard result for two-dimensional spacetime where the proper time is given 
by the relation r = (1 — 2M/r)^/^t. The remarkable thing about the above expression is the second 
term. It is known that there are three possible vacuums around a black hole that could be chosen 
for renormalization: the Boulware, Unruh and Hartle-Hawking vacuum states. However, the way 
in which the quantum inequality is developed picks out the Boulware vacuum. The mode functions 
that are used in deriving the quantum inequality are defined to have positive frequency with respect 
to the timelike Killing vector. This is also the condition for defining the Boulware vacuum state. 
The same is true for the four-dimensional black hole where the renormalized quantum inequality 
again selects the four-dimensional Boulware vacuum. 

Finally, in Chapter |^ we will look at the dynamic spacetime of the Alcubierre "warp drive" 
in the short sampling time limit. We will show for super luminal bubbles of macroscopically useful 
size that the bubble wall thickness must be on the order of only a few hundred Planck lengths |40|. 



Under this restriction, the total required negative energy (the energy density integrated over all 
space) is 100 billion times larger than the entire mass of the visible universe. This would seem to 
dash any hope of using such a spacetime for superluminal travel. 



It comes as no surprise that similar results are found for the Krasnikov spacetime |41, 42 1. 
In this spacetime, a starship traveling to a distant star at subluminal velocity, creates a tube of 
negative energy behind it, altering the spacetime as it goes. A good analogy here is the digging of 
a subway tunnel. When the ship reaches its destination it can then turn around and return home 
through the tube. To observers who remain stationary at the point of departure, the amount of 
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time that elapses between the departure and return of the spaceship can be made arbitrarily small, 
or even zero. Everett and Roman [42| showed that when the quantum inequality is applied to 
this spacetime in the short sampling time limit, the negative energy is constrained in a wall whose 
thickness is again only on the order of a few Planck lengths. 



Chapter 2 

How Negative Energies Arise in QFT 



In quantum field theory, there are certain states in which the local energy density is negative Q, 
thus violating the weak energy condition . One example is the vacuum plus two particle state for 
the quantized scalar field. Another is the squeezed states of light for the quantized electromagnetic 
field. We will look closely at both cases in Section |2.1| . 

Specific states in which the quantized fermion field can yield negative energy densities are 
just beginning to emerge. Vollick has recently demonstrated that the energy density can become 
negative when the particle content is a superposition of two single particle states traveling in 
different directions |^0|. In Section 2.2, we will demonstrate that the energy density for the fermion 
field is negative when the state is a superposition of the vacuum and a particle-antiparticle pair. 

In Section |2.3| , the negative vacuum energy density associated with the Casimir effect for a quan- 
tized field near perfectly refiecting boundaries will be discussed. For the quantized electromagnetic 
field, this gives rise to an attractive force between two parallel perfectly conducting planar plates. 
We will also look at how non-zero energy densities arise for the vacuum state in curved spacetimes. 



2.1 Quantized Scalar Field 

We begin our discussion of negative energy densities with the quantized scalar field in Minkowski 
spacetime. This model demonstrates many of the significant points about negative energy in 
quantum field theory. 

Let us assume that we have a massive, quantized scalar field (f){x) which resides in Minkowski 
spacetime and satisfies the Klein-Gordon equation 

- 5tV(x) + V V(x) - m2</>(x) = 0. (2.1) 

The field (j){x) can be Fourier expanded in terms of a complete set of plane wave mode functions 

h{x) = (2uL^y'^\''^-^-'^\ (2.2) 



each of which individually satisfies the wave equation (^^) for any choice of the mode label k = 
{kx, ky,kz). Note, we are using box normalization of the mode functions. The general solution to 
the wave equation above is given by 

(t>{x) = ^ [ak/k(x) + almx)] . (2.3) 



10 
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In the rest frame of an observer, the energy density is the time component of the stress-energy 
tensor, 



1 r 



m (j) [x 



(2.4) 



where we have chosen to use minimal coupling of the scalar field. Similar effects could be proven 
in general for any form of the coupling in Minkowski space. If we now insert the mode function 
expansion of the scalar field, we arrive at 



t 1 



\^{ + 1^ • ^0 /k/k'^k^k' - /k/k'^k^L ~ /k/k'Okak' + /k/k'^kflk^ 

kk' 

/k/k'^k^k' + /k/k'«k«k' + /k/k'ttkOk' + /k/k'^k^k' } 



+m 



(2.5) 



This particular equation has a divergence for any quantum state for which we may choose to 
evaluate its expectation value. The problem stems from the akoj^, terms of the above expression. 
For example, in the vacuum |0), the stress-tensor reduces to 



(0|r«|0) = \Y. H + k • k' + m2) 4/^, (0|akai,,|0) = ^ E (l^l' + 
kk' k 



m 



1/2 



(2.6) 



Here we have made use of the standard definitions of the creation and annihilation operators such 
that when they act on the vacuum state, they have the values 



4io) 



lu 



and ak|0) = 0, 



respectively. In addition, when they act on a number eigenstate they become 



Ok |"-k) = \/nTT |n + Ik) 



and 



Ok Ink) 



n \n 



(2.7) 



(2i 



We see once the summation is carried out in Eq. ( |2.6| ) that the answer is infinite. This divergence 
seems to indicate that the vacuum itself contains an infinite density of energy. In non-gravitational 
physics, absolute values of the energy are not measurable, only energy differences. It has become 
standard practice to renormalize the energy density (as well as the energy) by measuring all energies 
relative to the infinite background energy. Typically this is achieved by normal ordering such that 
all of the annihilation operators stand to the right of the creation operators within the expectation 
value. Then, we subtract away all of the divergent vacuum energy terms. For example, if we make 
use of the commutation relation for the akOk' terms in Eq. ( |2.5| ), we find 



kk' 



uJio' + k • k') 



/kA'Akflk' - /k/k'«k«k' ~ /k/k'Okak' + /k/k'^k'^^k 



+m 



/k/k'Ok^k' + /k/k'^k^k' + /k/k'OkOk' + /k/k'^k'^k] } + ^ 



The last term is the infinite vacuum energy, and is subtracted away to complete normal ordering. 
We define the normal ordered energy density as 



■■Tti 



Tti 



\Tu\0), 



(2.10) 



which produces a finite expectation value for physically allowable states. 
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2.1.1 Vacuum Plus Two Particle State 



We can now address how negative energy densities can be found. Let us consider the particle 
content to be defined by the superposition of two number eigenstates, e.g., 



f|0)+i|2.). 



If we look at the expectation value of the energy density we find, 

(^1 : Tu : IV 



UJ 

2l3 



\/6 

1 — cos 2(k • X — cot) 



(2.11) 



(2.12) 



where we have also assumed that the field is massless. We see that the energy density oscillates 
at twice the frequency of the mode that composes the state. We also see that once during every 
cycle the energy density becomes negative for a brief period of time. A plot of this is shown in 
Figure O. 




cot 



Figure 2.1: Time evolution of the energy density for the vacuum plus two particle state. Note 
that during every oscillation, the energy density dips below zero for a brief portion of the cycle. By 
manipulating the quantum state, it is possible to drive the energy density to increasingly negative 
values. Units on the vertical axis are in multiples of lo/2L^. 

Actually, there are an infinite set of vacuum plus two particle states that give rise to negative 
energy densities. To show this, let us begin with a general superposition of an n and an n + 2 state, 

IV') = c„|nk) + c„+2|n + 2k), (2.13) 

where the normalization of the state is expressed by the condition 

|c„|2 + |c„+2p = 1. (2.14) 

There is an additional restriction that \cn\ / 1. This corresponds to a pure number eigenstate 
because the unit normalization would cause Cn+2 = 0. 
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Figure 2.2: A plot showing the two conditions necessary to obtain negative energy densities for 
the vacuum plus two particle state. The normalization condition for the coefficients cq and C2 is the 
unit circle. The necessary condition for negative energy |co| > \/2|c2| is represented by the light 
gray region. The set of all possible states that meet both of these conditions is denoted by the bold 
curve at the top and bottom of the circle. The states with |co| = 1 must be excluded because they 
correspond to the vacuum state, where the energy density is zero. 
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For a superposition of the n and n + 2 particle states, the energy density is given by 



■■Ttt : 



I3 



n + 2|c„+2| 



UJ 



m 



[Re«c„+2) cos 2(k • X - ujt) - |Im(c;c„+2 



(n + l)(n + 2) X 

sin 2(k • X 



a; 



t)]. 



(2.15) 



In order to simphfy our calculations, let us take c„ and Cn+2 to be real. Then, in order to have 
negative energy densities, we must satisfy the condition 



> 



CO 



1 



uj^ - m? y/(n + l)(n + 2) 



2|c„+2| + 



n 



(2.16) 



In addition, we must satisfy the normalization condition ( p.l4| ), which for real coefficients is the 
equation for a unit circle in the (cn, Cn+2) plane. 

For the massless (m = 0) scalar field in the vacuum plus two particle state (n = 0) , the condition 
to obtain negative energy densities, Eq. ( p. 16 ), reduces to 



I Co 



> V2 



C2I 



(2.17) 



We can see from Figure 2.2 that any state, except |co| = 1, that lies on the bold sections of the circle 
will satisfy both the above condition to obtain negative energy densities and the unit normalization 
condition. For real coefficients, c„ and c„+2j only the vacuum plus two particle state can generate 
negative energy densities. If we plot the condition to find negative energy densities, Eq. ( 2.16| ), 
for Cn and 0^+2 real and n > 1, we find that the range of values that would allow negative energy 
densities is incompatible with the unit normalization condition. 



2.1.2 Squeezed States 

The squeezed states in quantum field theory can also produce negative energies. The squeezed 
states of light have been extensively investigated in the field of quantum optics and are realized 
experimentally. We will confine ourselves to the squeezed states of the quantized scalar field, 
although the treatment would be the same for the electromagnetic field. 



We begin with the definition of two additional operators. The first, introduced by Glauber |51| 
is the displacement operator 



Dk{z) = exp yza\^ 
which satisfies the commutation relations 



Z Ok 



4" ^U^) 



The second is the squeeze operator 

sue) = exp 

which satisfies the relations 

4(C)ak5k(C) = 



flk' + zS]^]i', 
al., + z*(5kk'- 



Z Ok 



ak cosh r — al^e^^ sinh r. 



aj^ cosh r 



a^e *'^sinhr. 



(2.18) 

(2.19) 
(2.20) 

(2.21) 



(2.22) 
(2.23) 
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The complex parameters z and ( may be written in terms of their magnitude and phases as 

z = se'^ and ( = re'\ (2.24) 

We can write a general squeezed state for a single mode as IQ 

\{z,Ch) = DUz)SUCm. (2.25) 

It has been shown that for a massless scalar field the expectation value of the energy density for 
the above squeezed state is given by 



P = {{z,C}k\-Ttt:\{z,Ch 
to 



^ |sinh r cosh r cos(26' + 5) + sinh^ r + [I - cos 2(6* + 7)] | , (2.26) 



where 9 = k-x— tjt. When C = 0; the squeezed state reduces to a coherent state and the expectation 
value of the energy density is always positive definite. The coherent states are interpreted as the 
quantum states that describe classical field excitations. Thus, the production of only positive energy 
in coherent states is consistent with classical sources for gravity obeying the WEC. 

A different case, known as the squeezed vacuum state, is defined when z = 0. Such states 
result from quantum mechanical particle creation. An example is second harmonic generation in 
non-linear optical media. For the squeezed vacuum state, the energy density is given by 

({0,C}k| : Ttt : |{0,C}k) = ^ sinhr [sinhr + coshr cos(2^ + 5)] . (2.27) 

The squeezed vacuum energy density will have the same negative energy density behavior as the 
vacuum plus two particle state above, with the energy density falling below zero once every cycle 
if the condition 

coshr > sinhr (2.28) 

is met. This happens to be true for every nonzero value of r, so the energy density becomes negative 
at some point in the cycle for a general squeezed vacuum state. In a more general state which is 
both squeezed and displaced, the energy density will become negative at some point in the cycle if 
s <C r. 

There are other interesting aspects of these states as well. When the quantum state is very 
close to a coherent state, we have seen that the energy density is usually not negative. In addition, 
fluctuations in the stress-tensor for these states are small compared to the expectation value of 
the stress-tensor in that state. On the other hand when the state is close to a squeezed vacuum 
state, there will almost always be some negative energy densities present, and the fluctuations in 
the expectation value of the stress-tensor start to become nearly as large as the expectation value 
itself. This seems to indicate that the semiclassical theory is beginning to break down 



2.2 The Quantized Fermion Field 

In this section we will show that it is possible to generate negative energy densities using the fermion 
fleld. Therefore it does appear possible to use fermions to violate the weak energy condition. One 
example of a fermion state that produces negative energies has been demonstrated by Vollick 
for a superposition of two single particle states traveling in different directions. The particular 
example that we will demonstrate has a particle content which is a superposition of the vacuum 
plus a particle-antiparticle pair. However the occurrence of negative energy densities for the above 
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particle states can be a frame-dependent phenomenon. If the particle-antiparticle pair is travehng 
in the same direction, but with differing momenta, then it is possible for the energy density to 
become negative. However, in the center of mass frame of the particle-antiparticle pair, the local 
energy density always has a positive definite value. 

The fermion field, denoted by ip{x), is the solution to the Dirac equation, 

h^V^tp{x) - mtpix) = 0, (2.29) 

where m is the mass of the spin-^ particle, and the 7^ are the the Dirac matrices in flat spacetime, 
which satisfy the anticommutation relations 

{7^ 7"} = 2^?'^^ (2.30) 
The general solution of ip{x) can then be expanded in plane wave modes as 

Here b{p, s) and its hermitian conjugate b^{p, s) are the annihilation and creation operators for the 
electron (fermion), respectively, while d{p, s) and its hermitian conjugate S{p, s) are the respective 
annihilation and creation operators for the positron. All four operators anticommute except in the 
case 

{b{p, s), b\p\ s')} = 6ss'6Hp - P') = {d{p, s), S{p\ s')] . (2.32) 

The annihilation and creation operators for spin-i particles satisfy the anticommutation relations 
as opposed to the commutation relations for the scalar field because the fermions must obey the 
Pauli exclusion principle. 

For the fermion field, the stress-tensor is given by [^3|| 

T 

~ 2 



^7(mV.)V' - (V(^V5) 7.)V'] • (2.33) 



Inserting the mode function expansion into the energy density, and then normal ordering by making 
use of the anticommutation relations, we find 



Ttt = . y / d?pd?p'^ ^ {E„ + E.A X 

2(27r)3 t^J' "^ 



\/EpEp 



bHp', s'Hp, s)u\p', s')u{p, s)e-^(P^-PM)^'' + S{p, s)d{p' , s')v\p', s')v{p, s)e+'^P^-P>' 
+^.Y.jd'pdy-^iE,-E,). 

X 'd{p', s')b{p, s)v\p', s')u{p, s)e-^(f^+PM)^'' - b^{p, s)S{p', s)u\p', s')v{p, s)e+'^P^+P''^>'' 
-J^Ejd'pE,. (2.34) 

Unlike the vacuum energy density for the scalar field, we see that the fermion field has an infinite 
negative energy density that must be removed by renormalization. Thus we again look at the fully 
normal ordered quantity, 

:Ttt:=Ttt- {0\Ttt\0). (2.35) 
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While we are about to show that there are specific states for which the energy density can become 
negative, one should note that the Hamiltonian, 



H = Jd'x: Tu ■■=J2j d^PEp [^Hp, s)b{p, s) + d^p, s)d{p, s) 



(2.36) 



is always a positive definite quantity. 

Now, we would like to show that the state 



\(f>) = coo|0) + c^tI^i *i ; ^2 S2) , 



(2.37) 



admits negative energies for some values of cqo and c^j. For purposes of clarity, we have defined 
the above states by 



\ks) = b\k,s)\0), 
= d^k,s)\0), 
\kiSi;k2S2) = d\k2,S2)b'^{ki,si)\0) = -\k2S2;kiSi) 



(2.38) 
(2.39) 
(2.40) 



If we take the expectation value of the energy density with respect to this vacuum plus particle- 
antiparticle state, we find 



r^m\Cii\ 



+ 



(2vr) 
m {El — E2 



2(27r)3 ^/WEi 



u^{ki,s-[)u{ki,si) + v\k2, S2)v{k2, S2) 

4oCii v\k2, S2Hki,si)e-'^'''+''^^^^ + C.C. 



(2.41) 



This expression is similar to the scalar field energy density for the vacuum plus two particle state. 
In order to have negative energy densities, we must find some combination of the mode functions 
and the coefficients cqo and such that the amplitude of the interference term in the expression 
above is larger than the magnitude of the first term, which is the sum of the energy of the particle 
and antiparticle. In order to proceed any further, we must choose a specific form for the mode 
functions. This involves choosing a representation of the Dirac matrices. Let, 



and 



.(...)^[M^.™),-''^(™;>), 



(2.42) 



(2.43) 



where p is the three momentum of the particle, a is the vector composed of the Pauli spin matrices, 
and (j){s) is a two-spinor. Using this basis, it is possible to show 



u' {k, s)u(k, s') = — Sgs' = vUk, s)v{k, s'). 
m 

In addition, the interference term between the particle-antiparticle states is 

0t(s2)<T • [ki{E2 + m)+ k2{Ei + m)] ^(si) 



v\k2,S2)u{ki,si) 



y/2m{Ei + m){E2 + m) 



(2.44) 



(2.45) 
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To evaluate the energy density explicitly at this point, let us take si = S2 = +1, therefore 

Hsi) = <Pis2)=(l). (2.46) 



^0, 

For simplicity let the propagation vectors ki and k2 both point in the z-direction, i.e., 

k'l = {Ei,0,0,ki) and k>^ = {E2,0,0, ki) , (2.47) 
and let cqo and c^i be real. Then the energy density can be written as 

^ ' (27r)3^ '^^2(27r)2 y^E^E^{E^+mj{E^+mj ^ ' " ^ ' 

By an appropriate choice of the momenta and the coefficients cqo and c^j this can be made negative. 
This is most easily seen in the ultrarelativistic limit when E ^ m. The energy density is 



(■ Tu :) = + .3 coocn cos(i?i + E2)t. (2.49) 



El + E2 2 El — E2 
(27r)3 + ~(2^ 

We then find that the energy density is negative, apart from the special case cqo = 1, c^j = 
corresponding to the vacuum, when the condition 

El + E2 



|coo| > 



-cii 



(2.50) 



El — E2 

is met. Thus, an infinite number of states for negative energy densities could possibly exist. 

It is interesting to note that in the center of mass frame of the particle-antiparticle pair, the 



local energy density in this state is a positive constant. This can easily be seen from Eq. (2.48|), or 



more generally from Eq. ( p.4lD , when Ei = E2 and ki = —k2. Thus, for a given particle content of 
the state, whether it is possible to detect negative energies is dependent upon the frame in which 
any measurement is to be carried out. 

2.3 The Vacuum Energy 

So far we have discussed only the cases where the particle content of the quantum state makes the 
energy density negative. These could be called coherence effects for the fields. As we have seen in 
the plane wave examples above, the negative energy densities are periodic in both space and time. 
For more general mode functions this need not be the case, and we could have various configurations 
of negative energy densities. However, it is also possible to generate negative energies without the 
presence of particles. The vacuum energy density of a spacetime can be positive or negative with 
respect to the Minkowski space vacuum after renormalization. Casimir originally proved this in 
1948 for the electromagnetic field between two perfectly conducting planar plates Q. Here the 
energy density between the two plates is negative. Similar effects can be shown in gravitational 
physics. For example, the vacuum energy is non-zero for the scalar field in Einstein's universe and 
can be positive or negative, depending on our choice of either minimal or conformal coupling for 
the interaction between gravity and the scalar field in the wave equation, 

9'"'y^Nv ^{x) - + iR{x)) 0(x) = 0. (2.51) 
Minimal coupling is ^ = 0, while conformal coupling is given by 
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where n is the dimensionahty of the spacetime. In four dimensions ^ = 1/6, while in two dimensions, 
minimal and conformal coupling are the same. The vacuum energy density for a massless scalar 
field in the four-dimensional static Einstein universe [15, 35] is 

0.411505 

Pvacuum — '. ^ a (Z.Oo) 



Pvacuum — + ^on^2„4 " (2-54) 



for minimal coupling. Meanwhile, for conformal coupling the vacuum energy density is ||13|, |14| 

1 

480T^a4 

Here, a is the scale factor of the closed universe. We will discuss more fully the derivation of the 
vacuum energy for minimal coupling in Einstein's universe in Section 4.3.3| , when we renormalize 
the quantum inequality in Einstein's universe. At this point we would like to discuss a slightly 
simpler example, the vacuum energy near a planar conductor in both two dimensions and four 
dimensions. 



2.3.1 Renormalization 

We begin with a massless scalar field in Minkowski spacetime with no boundaries. The positive 
frequency mode function solutions for the scalar wave equation are 

/k(x,t) = (4^a;)-i/2e^(k— (2.55) 

where — oo < k < +00 is the mode label, and to = |k| is the energy. In the case of arbitrary 
coupling, the stress-tensor is 

T^, = (1 - 2^0,^0,. + (2C - 1/2)7]^,^""^^ - 2e</',M.</' . (2.56) 

We could proceed by placing the mode function expansion for the field (;/)(x, t) into the expression for 
the stress-tensor above and then taking the expectation value as we did in the preceding sections. 
We would find the infinite positive energy density as we have previously and then subtract it away to 
find the renormalized energy density. There is an alternative method using the Wightman function, 
defined as 

D+{x,x') = {0\(l){x)(j){x')\0). (2.57) 

The expectation value of the stress-tensor above in the vacuum state can be written equivalently 
as 

(0|r^,|0) = hm [(1 - 20d^dl + (2^ - l/2)g^,gP^dpd'^ - 2^d^d,] D+{x,x'), (2.58) 

x'—yx 

where 9^ is the ordinary derivative in the unprimed coordinates and in the primed coordinates. 
The process of renormalization involves removing any infinities by subtracting away the equivalent 
infinity in the Minkowski vacuum. This can be accomplished by renormalizing the Wightman 
function before we act on it with the derivative operator to find the expectation value of the 
stress-tensor. The renormalized Wightman function is defined as 

Dtn.i^.o,') = D+{x,x') - I)J,,„fc_,,(x,x'). (2.59) 

We see that the renormalized vacuum energy in Minkowski space is then set to zero. A similar 
formalism can be carried out for curved spacetimes as well. To find the renormalized Wightman 
function, we would first find the regular Wightman function. Then we subtract away the Wightman 
function that would be found by taking its limit when the spacetime becomes fiat. However, in 
curved spacetimes, this does not remove all of the infinities in the stress-tensor. There may be 
logarithmic and/or curvature-dependent divergences which must also be removed. 



CHAPTER 2. HOW NEGATIVE ENERGIES ARISE IN QFT 



20 



2.3.2 Vacuum Energies for Mirrors 



Let us start by considering the case of a two-dimensional flat spacetime in which we place a perfectly 
reflecting mirror at the origin. In two-dimensional Minkowski spacetime, the Wightman function 
is found to be 



D 



Minkowksi 



{x, x') 



1 

47r 



In 



./\2 



it-t'Y-{z-z') 



(2.60) 



However, we do not expect this to be the Wightman function when the mirror is present. That 
is because the mode functions will be altered in such a way that they vanish on the surface of 
the mirror. Because the spacetime is still flat, just with a perfectly reflecting boundary, the new 
Wightman function can be found by the method of images, with a source at spacetime point to 
be x', and the observation point at x. The Wightman function with the mirror present is now 
composed of two terms, 



D+{x, x') = -i- In [(t - t'f -{z- z'f 



/^2 



{t - t'Y - (z + z') 



(2.61) 



To renormalize the Wightman function, we now subtract away the Minkowski Wightman function 
yielding 

(2.62) 



D 



1 



i?en.(^>^') = 4^1n [{t-t'f-{z + z'f 



The vacuum energy density in this spacetime is then 



(0 1 Tit 1 0) Ken. = lim 



\dtdt> + {\- 20d,d, 



47r 



{t - t'f -{z + z') 



l\2 



27rz^ 



2idtdt 

Similar calculations can be carried out for the other components of the stress-tensor to find 



. (2.63) 



(0|T^^|0)Ren. 



2t:z^ 



1 




(2.64) 



We see that the vacuum energy density is non-zero for any value of ^ 7^ 0. In addition, the energy 
density diverges as one approaches the mirror. The actual sign of the energy density depends on the 
sign of the coupling constant, so we can have infinite positive as well as negative energy densities. 
However, at large distances away from the mirror the vacuum energy density rapidly goes to zero. 

Similar results can be found for a planar perfectly reflecting mirror in four dimensions. With 
the mirror located in the xy-plane at z = the renormalized Green's function for a massless scalar 
field is 



Using Eq. ( |2.58|) we find 



1 

'4^ 



[x 



x'f + {y 



y) +{z + z^ 



f-{t- t'f 



(2.65) 



1-6^ 

{0\T^u\0)Ren. = ^g^2^4 
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10 




(2.66) 



As was the case with the two-dimensional plate, the four-dimensional vacuum energy still diverges as 
one approaches the mirror surface. In addition, there are different values of the coupling constant 
^ that cause different effects on the stress-tensor. For the coupling constant less than 1/6, the 
energy density is everywhere negative and the transverse pressures are positive. In addition, both 
are divergent on the mirror's surface. This includes the often used case of minimal coupling, ^ = 0. 
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When the coupling constant is equal to 1/6, the conformally coupled case, the energy density 
and transverse pressures vanish. Finally, when the coupling constant becomes greater than 1/6, 
the energy density becomes positive and is still divergent at the mirror. However, the transverse 



pressures now become negative. These are summarized in Table 2.1. One should note that for all 
values of the coupling constant, the component of the pressure that is perpendicular to the surface 
to the mirror, the z-direction, is always zero. 



Value of Coupling Constant 


Energy Density 


Transverse Pressure 


e 


(Ttt) 


{Txx) = {Tyy) 


e< 1/6 


< 


> 


^ = 1/6 








e> 1/6 


> 


< 



Table 2.1: Vacuum energy density and pressures for a perfectly reflecting mirror in four dimensions 
as a function of the coupling parameter ^. 



2.3.3 The Casimir Force 

Let us consider two perfectly conducting planar plates, parallel to each other and separated by a 
distance L along the z-axis. Between the two conducting plates, the stress-tensor of the scalar field 
for conformal coupling is given by 55 1 



(0|T^.|0) 



vr 



IJ,u\^/Conf. 



1440L4 

















1 














1 














-3 



(2.67) 



In the case of minimal coupling it is given by 



3 - 2sin2(7rz/L) 



(0|V|0)Mm. = {0\T^u\0)conf. + ■ 4, /. 

AaL^ sm (vrz/Lj 
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10 




(2.68) 



In the region outside of the plates, the stress-tensor is given by Eq. ( |2.66| ). The net force, F, 
per unit area. A, acting on the plates is then given by the difference of the Tzz component of the 
stress-tensor on opposite sides of the plate. For both minimal and conformal coupling, there is an 
attractive force between the two plates with a magnitude 



\FiL)/A\ 



vr 



480 ' 



(2.69) 



Similarly, for the quantized electromagnetic field between two conducting plates, the attractive 
force between the plates is 



\F{L)/A\ 



(2.70) 



240 ■ 

The factor of two difference comes about because the electromagnetic field has two polarization 
states. Lamoreaux has confirmed the existence of this force experimentally |1Q]. This was done by 
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measuring the force of attraction between a planar disk and a spherical lens, both of which were 
plated in gold. The Casimir force is known to be independent of the molecular structure of the 
conductors, provided the material is a near perfect conductor, but sensitive to the actual geometry 
of the plates. Attempts to measure the force between two parallel plates were unsuccessful because 
of the difficulty in trying to maintain parallelism to the required accuracy of 10"^ radians. For the 
plate and sphere, there is no issue of parallelism, and the system is described by the separation of 
the points of closest approach. However the Casimir force for this geometry is not given by (|2.70| ), 
but has the form 



Fia) 



27raR 



3 y240a4 J ' 



(2.71) 



where R is the radius of curvature of the spherical surface and a is the distance of closest approach. 
In addition, the force is independent of the plate area. Lamoreaux finds that the experimentally 
measured force agrees with the predicted theory at the level of 5%. With direct confirmation 
of the Casimir force, we must admit that the vacuum effects in the stress-tensor are not merely a 
mathematical curiosity, but have direct physical consequences. In general relativity this is extremely 
important because the vacuum stress-tensor can act as a source of gravity. 



2.3.4 Vacuum Stress-Tensor for Two-Dimensional Stars 

In all of the examples of the preceding section for vacuum energy densities, we studied how the 
effects of boundary conditions on the scalar field can induce the vacuum to take on a non-zero 
value in flat spacetime. We saw that the energy density near a perfectly reflecting mirror becomes 
divergent on the surface to the plate. It is interesting to note that non-zero vacuum energies can 
also be found in gravitational physics. In some sense, the non-zero value of the stress-tensor is more 
important in gravitational physics because it can serve as a source in Einstein's equation. Thus 
the vacuum energy can have non-trivial effects on the evolution of the spacetime. The most well 
known, and well studied example of this is the evaporation of black holes, originally demonstrated 



by Hawking |5£]. In this section we will study a slightly different problem, the vacuum energy for 
a massless scalar field on the background spacetime of a constant density star. 

We begin with the metric of a static spherically symmetric star composed of a classical incom- 



pressible constant density fluid. The metric for this type of star is well known |57| 



-B{r)df + A{r)dr'^ + r'^d^'^ 



(2.72) 



where dO^ is the length element on a unit sphere, and the functions A{r) and B{r) are defined by 



A{r) 



[1 - 2MrV^^ 
[1 - 2M/r]"^ 



for r < R, 
for r > R, 



(2.73) 



and 



B(r) 



(1 - 2Mr^/R 



3^1/2 



for r<R, 
for r > R. 



(2.74) 



\ 3(1 - 2A//i?)^/^ 
[1 - 2M/r] 

Here M is the mass of the star and R its radius. This metric has one constraint: The mass and 
radius must satisfy 

R>- (2M) . (2.75) 
8 

Otherwise the pressure at the core of the star will become singular. This radius is just slightly 
larger than the Schwarzschild radius of a black hole of the same mass. When a star becomes 
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sufficiently compact and its radius approaches the above limit, then the gravitational back-reaction 
should start to have a significant effect on the spacetime. For the rest of our calculations we will 
presume that condition ( 2.75 ) is satisfied. It would be a rather arduous task to find the vacuum 
stress-tensor for the four-dimensional spacetime. However, we can make significant progress for the 
two-dimensional equivalent with the metric 



ds^ = -B{r)dt^ + A{r)dr^ 
This metric can be recast into its null form 

ds^ = -B{r)dUdV, 
where the null coordinates are defined by 

U = t-r* and V = t 

with r* given by 



lAjr) 
B(r) 



dr. 



(2.76) 



(2.77) 



(2.78) 



(2.79) 



It is clear from the form of the metric (2.77) that this spacetime, like all two-dimensional spacetimes, 
is conformally fiat. Thus, we can apply the method of finding the vacuum stress-tensor in two- 
dimensional spacetimes, [ ^ ] and Section 8.2 of |53]. The stress-tensor components of the vacuum, 
obtained when positive frequency modes are defined by the timelike Killing vector, in the (U, V) 
coordinates are 



and 



{0\Tuu\0) = {0\Tvv\0) 



{0\Tuv\0) 



1 



1927r 
1 



2BB" - (F 



/\2 



967r 



BB" 



(2.80) 



(2.81) 



where the prime denotes the derivative with respect to r. By a coordinate transformation, we can 
calculate the components of the stress-tensor in the original (r, t) coordinates, 



(0|r^,(t,r)|0) 



1 

24^ 



B{r)B"{r) - \B\rf 











\A{r)B'{rf/B{r) 



(2.82) 



We can now directly calculate the value of the vacuum stress-tensor both outside and inside the 
star. In the exterior region of the star, the stress-tensor is given by 



(0| -^itj' I O)ea;terior — r\A 

z47r 











-(MVr^)(l - 2M/r) 



(2.83) 



This is the well known result for the Boulware vacuum of a black hole . We can see that for the 
stationary observer at a fixed radius outside of the star that the energy density is 



Pexterior 

(r) 



(7MVr'' - 4M/r3) 
247r(l - 2M/r) 



(2.84) 



For the case of the star, the exterior region does not reach all the way to the Schwarzschild radius, 
so there is no divergence as there was for a black hole. However, the energy density for an observer 
fixed at some constant distance away from the star is still negative. 
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The situation is somewhat different in the interior of the star. The stress-tensor components 
are given by 



(0|V|0) 



fiu I ^ / interior 



A{r) 
247r 



B{r)§ 



3^^(0 (l - - 1 + 






-A{r)Mlrl 



(2.85) 



The vacuum energy density for the interior of the star as seen by a stationary observer is 

-1 



PinteriorKf ) - I i ^3 



M\ 

R^ 



3iM(r) 1 



2M\ Mr2 



(2.86) 



It is not immediately obvious, but the energy density in the interior is everywhere positive, and 



grows in magnitude as we move outward from the center of the star, as plotted in Figure |2.3| . In 
contrast, the vacuum energy outside the star is negative and identical to that of the black hole in 
the Boulware vacuum state. The discontinuity in the vacuum energy is an effect of the star having 
a very abrupt change in the mass density when crossing the surface of the star. We presume that if 
the star makes a continuous transition from the constant density core to the vacuum exterior then 
the vacuum energy would smoothly pass through zero. 



PvacW 



0.00075 
0.00050 
0.00025 


-0.00025 
-0.00050 
-0.00075 



Figure 2.3: The vacuum energy density induced inside and outside a star. Here we have taken the 
star to be composed of some form of classical matter that has a constant density throughout. The 
vacuum energy is then calculated for a scalar field which does not interact with the matter, except 
that it is subject to the spacetime curvature. Units have been normalized such that the mass of 
the star is unity. The radius of the star is 4 units. 

The total vacuum energy of this spacetime can be found by integrating the local energy density 
over all space. The total positive vacuum energy contained inside the star is 



rR II 2M \ 

Einside = 2 p{r) dr = -^-^ lOM - V2MR tmilC^ V ' (^'^^^ 
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The contribution outside the star gives a total finite negative energy 

Eoutside = 2 / p{r)dr = --— 
JR 96tt 

The total energy due to the vacuum is the sum 



14M 2 1 / 2M\ 
In 1 I 



ii2 R M \ R J 



(2i 



Erotal — Einside + Eoutside ■ (2.89) 



This is always finite and negative, as shown in Figure 2.4 along with the individual contributions 
from the interior and exterior of the star. For very diffuse stars, R ^ 2M, the total positive interior 
vacuum energy very nearly compensates for the negative energy outside the star. In this limit the 
leading contributions to the total energy go as 

i^To,.~-:^-5^-0(fl-)--. (2.90) 

In the other regime, when the star becomes very compact, the Boulware vacuum begins to dominate 
and in the limit R — > 9M/4, the energy density goes to its finite minimum value -Emm = —5. 306574 x 
10~^M^^ This residual energy can be considered a mass renormalization of the star due to the 
vacuum polarization, 

5 306574 X lO"'^ 

M = Mbare + 6M = Mbare " ^PZancfc- (2-91) 



Unless the mass of the star is very close to the Planck mass, Mpianck, this correction to the effective 
mass of the star is extremely small, accounting for at most half a percent of the star's mass. 
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-0.005 



-0.010 




Figure 2.4: The total integrated vacuum energy density induced inside and outside a star as a 
function of the radius. The top line is the contribution from the positive energy in the interior of 
the star. The lower line is the contribution from the exterior Boulware vacuum energy. The middle 
line is the total over all space. The energy is in units of inverse mass. 



Chapter 3 

Derivation of the Quantum 
Inequalities 

3.1 General Theory for Quantized Scalar Fields 

We begin by considering the semiclassical theory of gravity, where the classical Einstein tensor 
on the left-hand side of Eq. ( |3.1| ) is equal to the expectation value of the stress-energy tensor 
(stress-tensor) of a quantized field on the right-hand side, 

G,,, = 8^ (V) . (3.1) 

Here we are using Planck units, in which h = c = G = 1. We will take the stress-tensor to be that 
of the massive, minimally coupled scalar field, (j){x), given by 

Tf,u = (t>;fi(t>;iy " ^fl'M'^'A'^'A;^ + ^"l^^^i^c/i^ , (3.2) 

where (p-a denotes the covariant derivative of (j) on the classical background metric and m is the 
mass of the field. We shall develop quantum inequalities in globally static spacetimes, those in 
which dt is a timelike Killing vector. Such a metric can be written in the form 

ds'^ = -|5«(x)|dt2 + gij{^)dx'dx^ , (3.3) 

where the function 5ti(x) is related to the red or blue shift dependent only on the observer's position 
in space, and gij{x) is the metric of the spacelike hypersurfaces that are orthogonal to the Killing 
vector in the time direction. With this metric, the wave equation 

□0-m20 = ^= fa«.A^5"^a/30) -^20 = (3.4) 



becomes 

- T^df(l) + V*Vi0 - = 0, (3.5) 
\9tt\ 

where g = det{g^u). The timelike Killing vector allows us to use a separation of variables to find 
solutions of the wave equation. The positive frequency mode function solutions can be written as 

/A(x,t)= Ux{^)e~''^\ (3.6) 
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The label A represents the set of quantum numbers necessary to specify the mode. Additionally, 
the mode functions should have unit Klein-Gordon norm 



{fx, fx') = ^XX'- 



(3.7) 



A general solution of the scalar field (p can then be expanded in terms of creation and annihilation 
operators as 

</' = E(«A/A + al/;), (3.8) 

A 

when quantization is carried out over a finite box or universe. If the spacetime is itself infinite, 
then we replace the summation by an integral over all of the possible modes. The creation and 
annihilation operators satisfy the usual commutation relations []53| . 

In principle, quantum inequalities can be found for any geodesic observer [^]. In many curved 
spacetimes, static observers often view the universe as having certain symmetries. By making use 
of these symmetries, we can often simplify calculations or equations. To observers moving along 
timelike geodesies the symmetries of the rest observers may not be "observable." In the most general 
sense, the mode functions of the wave equation in a moving frame may become quite complicated. 



Thus, we will concern ourselves only with static observers, whose four- velocity, = (- 



'\9tt\ 



,0), is 



in the direction of the timelike Killing vector. The energy density that such an observer measures 
is given by 



1 



\gtt\ 



1 



\9tt 



2 ±2 

m (p 



(3.9) 



Upon substitution of the mode function expansion into Eq. (p.9|), we find 
P = Re^ 



AA' 
+ 



LVLV 
\9tt\ 



lay+V^UxVjUye 



-i{u!+Lj')t 



axay 



[u*xUye 



a[ay + UxUye 



-i{uj+Lo')t 



axay 



+ 2 E ( ^^^^f^^ + ^'UlV.Ux + m^UlUx 



(3.10) 



The last term is the expectation value in the vacuum state, defined by a^lO) = for all A, and 
is formally divergent. The vacuum energy density may be defined by a suitable regularization and 
renormalization procedure, discussed in more detail later in Sections 4.2.3| and 4.3.2. In general, 
however, it is not uniquely defined. This ambiguity may be side-stepped by concentrating attention 
upon the difference between the energy density in an arbitrary state and that in the vacuum state, 
as was done by Ford and Roman |^2|. We will therefore concern ourselves primarily with the 
difference defined by 

:p:=p-(0|/5|0), (3.11) 

where |0) represents the Fock vacuum state defined by the global timelike Killing vector. 

The renormalized energy density as defined above is valid along the entire worldline of the 
the observer. However, let us suppose that the energy density is only sampled along some finite 
interval of the geodesic. This may be accomplished by means of a weighting function which has a 
characteristic time to- The Lorentzian function. 



(3.12) 
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is a good choice. The integral over all time of f{t) yields unity and the width of the Lorentzian is 
characterized by to. Using such a weighting function, we find that the averaged energy difference 
is given by 



TT 



{: Tu/\gtt\ ■.)dt 
t' + tl 



XX' 



\gtt 



+ m2 [c/^*f/ve-l--'l*°(alay) + C/AC/A'e-(-+-')*«(aAav)]} . 
We are seeking a lower bound on this quantity. It has been shown that 
Re^/(A)*/(A')e-l'^-"'l*°(aiav) > ReE/W*/(A')e-("+"'^*"(4«A'>- 

AA' AA' 

Upon substitution of this into Eq. ( 3.13| ) we have 



> Re^ 

AA' 



LJLJ 



\9tt\ 



UxU\'{a\ax') - UxUyiaxay) 



+ [v^U*xVjUy{a{ax') + V^UxV jUy{axay 



+m 



UlUy{a{ay) + UxUyiaxay)] } e-("+'^')*o . 



(3.13) 



(3.14) 



(3.15) 



We may now apply the inequalities proven in Appendix^. For the first and third term of Eq. (|3.15| ), 



apply Eq. ( |A.7| ) with hx = / ^/\gtt \ Uxe and hx = mUxe respectively. For the second term 
of Eq. ( p^ , apply Eq. (^) with Aij = gij and h\ = V'Uxe''^^'' ■ The result is 



\gtt\ 



U*xUx + V^U*x V,Ux + m'U*xUx e 



-2aJAio 



(3.16) 



This inequality may be rewritten using the Helmholtz equation satisfied by the spatial mode func- 
tions: 

V^VjUx + {u^Vlattl - m^)Ux = , (3.17) 

to obtain 



> - ^ ( ^ + iv^V,- ) iC/Ape-^-Ato 



This can be rewritten as 



Ap > 



1 / d] 



to 



4 V \9tt\ 



+ V^vAJ2\Ux{^)\'e 



(3.18) 



(3.19) 



There is a more compact notation in which Eq. ( |3.19| ) may be expressed. If we take the original 
metric, Eq. (|3.3|) , and Euclideanize the time by allowing t ^ Uq, then the Euclidean box operator 
is defined by 

d?. 



-'to 
\9tt\ 



+ V^V 



J • 



(3.20) 
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The analytic continuation of the Feynman Green's Function to imaginary time yields the Euclidean 
two-point function. The two are related by 

G£;(x, t; x', t') = iGpi^, —it; x', —it'). (3-21) 

In terms of the mode function expansion, the Euclidean Green's function is given by 

Gi,(x,-to;x,+to) = Elf^A(x)|V2-^*°, (3.22) 

A 

where the spatial separation is allowed to go to zero but the time separation is 2tQ. 

This allows us to write the quantum inequality for a static observer in any static curved space- 
time as 

A/5> -^□i5Gij(x,-to;x,+to). (3.23) 

Given a metric which admits a global timelike Killing vector, we can immediately calculate the 
limitations on the negative energy densities by either of the two methods. If we know the solutions to 
the wave equation, then we may construct the inequality from the summation of the mode functions. 
More elegantly, if the Feynman two-point function is known in the spacetime, we may immediately 
calculate the inequality by first Euclideanizing and then taking the appropriate derivatives. 

It is important to note that while the local energy density may be more negative in a given 
quantum state than in the vacuum, the total energy difference integrated over all space is always 
non-negative. This follows because the normal-ordered Hamiltonian, 

■.H := f :Tu: V^rx = Y.uJxa{ax, (3.24) 

A 

is a positive-definite operator; so {: H :) > 0. 



3.2 Quantum Averaged Weak Energy Condition 



Let us return to the form of the quantum inequality given by Eq. (3.18), 



^P^-T.(j^\+ i^J- ) |t^A(x)pe-2--*« . (3.25) 



A 



Since we are working in static spacetimes, the vacuum energy does not evolve with time, so we can 
rewrite this equation simply by adding the renormalized vacuum energy density Pvacuum to both 
sides. This is the vacuum in which the mode functions are defined to have positive frequency with 
respect to the timelike Killing vector. We then have 

PRen. > - E ( + I^'^i 1 |f^A(x)|2e-2-^*» + /9„™(x) , (3.26) 



A 



\9tt\ 



where pnen. is the sampled, renormalized energy density in any quantum state. Taking the limit of 
the sampling time to ~^ oo, we find (under the assumption that there exist no modes which have 
ujx = 0) that 

to f°° {Ttt/\gtt\)Ren.,, ^ , N o'7\ 

hm — / o at > pvacuumi^)- (3.27) 

to-*oo vr J-oo t + to 
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This leads directly to the Quantum Averaged Weak Energy Condition for static observers |32| 

f+OO 



{{lp\Ttt/\gu\\tp)Ren. - Pvacuum) dt > 0. (3.28) 

This is a departure from the classical averaged weak energy condition, 

/ + 00 
{i^\Tu/\gtt\ \i')Rendt > 0. (3.29) 
'OO 

We see that the derivation of the QAWEC leads to the measured energy density along the observer's 
worldline being bounded below by the vacuum energy. 

Recently, there has been considerable discussion about the vacuum energy and to what extent 
it violates the classical energy conditions. For example, Visser looked at the specific case of the 
violation of classical energy conditions for the Boulware, Hartle-Hawking, and Unruh vacuum states 



|18, 1£, 20, 21] around a black hole. However the vacuum energy is not a classical phenomenon, so 
we should not expect it to obey the classical energy constraints. From the QAWEC we see that the 
sampled energy density is bounded below by the vacuum energy in the long sampling time limit. 



3.3 Expansion of the QI for Short SampUng Times 



We now consider the expansion of the two-point function for small times. We assume that the 
two-point function has the Hadamard form 



G{x,x) - ^ 



AV2 



+ V\n{a + ie) + W 



(3.30) 



a + le 

where 2o"(x,x') is the square of the geodesic distance between the spacetime points x and x', 

A = -g-^'\x) dei{a,ab') g-^'\x), (3.31) 

is the Van Vleck-Morette determinant, and V{x^x') and W{x,x') are regular biscalar functions. In 
general, these functions can be expanded in a Taylor series in powers of o" |60|, 



V{x,x') = ^K(x,x')(T^ 

n=0 



where V^, and similarly are also regular biscalar functions with 

Vq = Wo - ^^aO-" + ^t^OafeO-'^CT^ + ^(-^VOafe;c + J^^0;(afec))cr''o-'''7'' + 



Vl 



2 ' 



and = a''^ . The coefficients, vq^ voab, ■ ■ ■ are strictly geometrical objects given by 



^^0 



VOab 



(3.32) 

(3.33) 
(3.34) 

(3.35) 
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1 



)R + 
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1 . 
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-)ni? + -(e 



-m 



(3.36) 



(3.37) 



We can then express the Green's function as 

1 + i^iJab^V^ 



G(x, x) 




a + ie 
+(^1 - 2^i;«^ + 



+ 



+ ...) 



ln(fT + ie) + 



(3.38) 



where we have also used the Taylor series expansion of the Van Vleck-Morette determinant 

= 1 + - ^Rah-rCJ^a^'cj'' + • • • . (3.39) 



The state-dependent part of the Green's function, W , is neglected because it is regular as fi ^ 0. 
The dominant contributions to the quantum inequality come from the divergent portions of the 
Green's function in the o" — > limit. 

We must now determine the geodesic distance between two spacetime points, along a curve 
starting at (xq, —to) and ending at (xq, +to)- For spacetimes in which \gtt\ = 1, the geodesic path 
between them is a straight line. Therefore, the geodesic distance is simply 21^. However, in a more 
generic static spacetime where gtt{^) is not constant, the geodesic path between the points is a 
curve, with the observer's spatial position changing throughout time. Thus, we must now solve 
the equations of motion to find the geodesic distance between the spacetime points. In terms of an 
affine parameter A, the geodesic equations are found to be 

dt at 



and 



ff«(x(A))| 

(fx^ : dxf^ dx" 
~d)^^ dX dX 







(3.40) 



0. 



(3.41) 



where at is an unspecified constant of integration. The Christoffel coefficients are 

1 



{dmj.k ~l~ 9mk,j 9jk,m) • 

It is possible to eliminate A from the position equations, and write 



^9'^ \9u\j 
0, 

2^ 



d'^x' 



1, 



dt^ +2'^"' 



dx^ dx^ 
+ T\k h 

^'^ dt dt 



\gtt\,k dx' dx^ 
\gtt\ dt dt 



0. 



(3.42) 



(3.43) 



Now if we make the assumption that the velocity of the observer moving along this geodesic is 
small, then to lowest order the second term can be considered nearly constant, and all the velocity- 
dependent terms are neglected. It is then possible to integrate the equation exactly, subject to the 
above endpoint conditions, to find 



x\t) 



^ IS'tt lx=Xo 



it' 



tn) + x: 



■ 



(3.44) 



CHAPTER 3. DERIVATION OF THE QUANTUM INEQUALITIES 



33 



We see that the geodesies are approximated by parabolae, as would be expected in the Newtonian 
limit. A comparison of the exact solution to the geodesic equations and the approximation is shown 
in Figure 3J- for the specific case of de Sitter spacetime. We see that the approximate path very 
nearly fits the exact path in the range —to to +to- 

The geodesic distance between two spacetime points, where the starting and ending spatial 
positions are the same, is given by 



As 



-to 



to 



\9tt (x(t))| +g,,(x(t))— — dt. 



(3.45) 



In order to carry out the integration, let us define 



f(t) = J-\gu (x(t))|+ff,j(x(t)) 



dt dt 



(3.46) 



We can expand f{t) in powers of t centered around t = 0, and then carry out the integration to 
find the geodesic distance. The parameter a can now be written as 



<7(xo,to) 



cr(xo, -to; Xo, +to) = 2 



V(o)/(^^)(o) + i/"(o)2 

5 3 



(3.47) 



However, we do not know the values of the metric at the time t = 0, but we do at the initial or final 
positions. Therefore, we expand the functions /(O) around the time —Iq. Upon using Eq. ( 3.44| ), 
we find 

1 ... .... 



<T(xo,to) ~ -2|5ti(xo)| ^0 - 7fi'tt'*(xo)fi'tt,i(xo)to + 



and 



o-*(xo,to) ~ 2to + 



1 gtt''{y^o) 9tt,i{y^o) 



4 + 



(3.49) 



3 |5«(xo)| 

In any further calculations, we will drop the xq notation, with the understanding that all of the 
further metric elements are evaluated at the starting point of the geodesic. Using Eq. (3.21) we can 
then write the Euclidean Green's function needed to derive the quantum inequality in increasing 
powers of Iq as 



GE{x,to) = ^ 



1 - 0{tl) + 



+ vq ln(2|5ftt| tl - ^gtt''gtt,i 4) 



4 - latv'att.i tl 

+ {vQ,k \gtt[^ + '^vi\gu\ - 2vooo) tl\D.{2\gu\ 4 - ^gtt'^gtt,i 4) ^ 



(3.50) 



Notice that none of the geometric terms, such as vq, changes during Euclideanization because they 

Ap > ( -^dl + ) Gij(x, -to; x, +to). (3.51) 



are time independent. The quantum inequality, ( 3.23 ), can be written as 



\9tt\ 
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Figure 3.1: The exact geodesic path (dotted hne) and the parabohc approximation (sohd hne) in 
de Sitter spacetime. The coordinate distance from r = to the horizon is a in the static coordinates. 
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If we insert the Taylor series expansion for the EucUdean Green's function into the above expression 



and collect terms in powers of the proper sampling time tq, related to to by tq 
write the above expression as 



\9tt\ 



0, we can 



Ap > 



327r2' 







1 + 



1 



m 



4 s^.'- 

3 V6 9tt 



Idttl 



T^H2r^) + 0{4) + 



In the limit of tq 



0, the dominant term of the above expression reduces to 

3 



A/5 > 



327r2ro^ 



(3.52) 



(3.53) 



which is the quantum inequality in four-dimensional Minkowski space p2| , Thus, the term in 
the square brackets in Eq. ( 3.52| ) is the short sampling time expansion of the "scale" function |35|, 
and does indeed reduce to one in the limit of the sampling time tending to zero. The range of 
sampling times for which a curved spacetime can be considered "roughly" flat is determined by the 

If this happens to be the Tq term, then for the first 



first non-zero correction term in Eq. ( |3.5 
correction to be small compared to one implies 



To < 



-1 



9tt V*Vi 



m 



-1/2 



(3.54) 



Each of the three terms on the right-hand side of this relation has a different significance. The m? 
term simply reflects that for a massive scalar field, Eq. ( p. 52 ) is valid only when the sampling time 
is small compared to the Compton time. If we are interested in the massless scalar field, this term is 
absent. The scalar curvature term, if it is dominant, indicates that the flat space inequality is valid 
only on scales which are small compared to the local radius of curvature. This was argued on the 
basis of the equivalence principle by various authors |44, 40, but is now given a more rigorous 



demonstration. The most mysterious term in Eq. (3.54) is that involving gtt- Typically, this term 
dominates when the observer sits at rest near a spacetime horizon. In this case, the horizon acts as 



a boundary, so Eq. ( 3.54 ) requires tq to be small compared to the proper distance to the boundary. 
In the particular case of \gtt\ = 1, we have tq = to and Eq. ( ^.52 ) reduces to 



Ap > 



327r2t4 



l+'-C-R 
3 V6 



m 



1 

+ 3 



12 



+ 4^;i - 4^;ooo ioM2to) + 



(3.55) 



This result has also been obtained by Song [31|, who uses a heat kernel expansion of the Green's 
function to develop a short sampling time expansion. We can now apply Eq. ( 3.55| ) to a massless 
scalar field in the four-dimensional static Einstein universe. The metric is given by 



-dt' 



dx^ + sin^ xide"^ + sin^ Od^"^) 



and the scalar curvature i? = 6/a^ is a constant. It can be shown that vi 
leads to a quantum inequality in Einstein's universe of the form 



^^000 



Ap > 



(3.56) 

l/8a^. This 



(3.57) 



In Section [4.3.2| , an exact quantum inequality valid for all to /a will be derived. In the limit to <C a, 
this inequality agrees with Eq. ( 3.57| ). Similarly, the exact inequality for the static, open Robertson- 



Walker universe will be obtained in Section 4.3.1, and in the limit to <^ a agrees with Eq. (3.55) 
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3.4 Electromagnetic Field Quantum Inequality 

In this section, we derive a quantum inequality for the quantized electromagnetic field in a static 



spacetime. It has been shown |62, ^ that the covariant Maxwell's equations for the electromagnetic 
field in a curved spacetime, 

Ff".^ = (3.58) 

and 

Ff,u;a + F^„,^ + F„^,^ = , (3.59) 

can be recast into the form of Maxwell's equations inside an anisotropic material medium in Carte- 
sian coordinates, by using the constitutive relations 

A = eifc^fc-(GxH)i (3.60) 

and 

Bi = iJikHk + (G X E)i , (3.61) 

where 

eik = IJ'ik = - {-gf^— and Gi = - — . (3.62) 
500 500 

Here the effects of the gravitational field are described by an anisotropic dielectric and permeable 
medium. However, when we consider the metric 

ds"^ = -\gtt{^)\dt^ + gij{^)dx^dx^ , (3.63) 

there is considerable simplification. Because g^i = for all i, the vector Gi is always zero. The 
constitutive relations are then simply given by 

D = eE and B = eH, (3.64) 

where 





B = 


eH 


11 


g'' 


^13 


21 




^23 


31 


g32 


^33 



6 = 6(X) = ^ I ^23 I . (3.65) 

ml 

The source-free Maxwell equations, in terms of (E,H), are given by 

VxE = -^, V-D = 0, (3.66) 

VxH = ^, V-B = 0. (3.67) 

We may also define the source-free vector potential yl^(x, t) = (0,A) with the relations to the 
electric and magnetic fields given by 

dA 

E = — — and B = V x A. (3.68) 

It is straightforward to show that the vector potential satisfies the wave equation 

52 A 

e^^V X r\V X A) = . (3.69) 
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The left-hand side of Eq. ( ^.691 ) involves only derivatives with respect to the position coordinates, 
and the right-hand side, temporal derivatives. This is a clear separation of variables, allowing us 
to write the positive frequency solutions as 



fk\x,i) = U^(x) e 



-iuit 



(3.70) 



where k is the mode label for the propagation vector and A is the polarization state. The vector 



functions, U^(x), are the solutions of 



r^V X e-^(V X U^) - w^u^ = 0, 



(3.71) 



and carry all the information about the curvature of the spacetime. The mode functions for the 
vector potential are normalized such that 



Ik 1 ^k' 



i^kk'^AA' • 



The general solution to the vector potential can then be expanded as 

A(x,t) =^(akAU^(x)e— * + aLur(x)e+-*) . 



(3.72) 



(3.73) 



k,A 



When we go to second quantization, the coefficients, al^^^ and OkA become the creation and annihila- 
tion operators for the photon. We will again develop the quantum inequality for the electromagnetic 
field for a static observer, as was done above for the scalar field. The observed energy density, p, is 
given by @ 

P=li-g)-'/HB-n + B.U). (3.74) 



Upon substitution of the mode function expansion, and making use of constitutive relations ( 3.64 ) 
we find 



1 



P 



Re J2 [4A«k'A' (U^I 9'' U^',) e*^^-"')* - OkAOk'A' (U^, g'' V^'j) e 



■i(LLl+U)')t 



1^**1 kk'.AA' 

+\gu\ Re [4A«k'A' (V X g^^ (V x U^',), e*(--')* 

kk'.AA' 

+akAak'A' (V X U^), g'^ (V x U^',), e-^(-+-')*" 
f^(UL g'^ + \gu\ (V X U^), g'^ (V x uf), 



1 V- 



k,A 



(3.75) 



The last term of the above expression is the vacuum self-energy of the photons. As was the case 
for the scalar field, we will look at the difference between the energy in an arbitrary state and the 
vacuum energy, i.e., 

:p:=p-(0|HO). (3.76) 
Again we integrate the energy density along the worldline of the observer, weighted by the Lorentz- 



ian sampling function. We may then apply the first inequality, Eq. (3.14) proven in |^8[. To the 
result we apply the inequality proven in Appendix We then find that the difference inequality 
on the energy density for a quantized electromagnetic field is given by 



k,A 



\gtt 



(uL g'' u^*) + b"l (V X U^), g'= (V X 



(3.77) 
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This expression is similar in form to the mode function expansion of the scalar field quantum 
inequality, and also reduces to an averaged weak energy type integral in the infinite sampling time 
limit. As was the case for the scalar field, the electromagnetic field quantum inequality ( 3.77D tells 



us how much negative energy an observer may measure with respect to the vacuum energy of the 
electromagnetic field. In order to find the absolute lower bound on the negative energy density we 
would have to add the renormalized vacuum energy into the above expression. 

This quantum inequality can be easily evaluated in Minkowski spacetime. In a box of volume 

with periodic boundary conditions, the mode functions are given by 

Vi, = (lujL^y^'^ e^^-"" e^, (3.78) 



where is a unit polarization vector and uj = Vk • k. Inserting the mode functions into Eq. ( 3.77] ) , 
and using the fact that k • = 0, we find 

^^^-^E^e-2-*o. (3.79) 

k,A 

The summation over the spin degrees of freedom yields 

Ap> --3_^^e-2^*«. (3.80) 
^ k 

In the continuum limit, L oo, the vacuum energy density vanishes, and the renormalized quantum 
inequality is found to be 



1 



oo 



oo 



167r2 t^- 



(3.81) 



This quantum inequality for Minkowski spacetime was originally proven by Ford and Roman 
using an alternative method. Comparison with the quantum inequality for the scalar field in 
Minkowski space, Eq. (|1.13| ), shows that the electromagnetic field quantum inequality differs by a 
factor of 2. This is a result of the electromagnetic field having two polarization degrees of freedom, 
unlike the scalar field which has only one. 



Chapter 4 

Scalar Field Examples 



4.1 Two-Dimensional Spacetimes 

There are a number of interesting results for two-dimensional spacetimes that we will discuss. The 
unique conformal properties in two dimensions allow the quantum inequality for all two-dimensional 
static spacetimes to be written in the form 



1 

^0 



where tq is the proper time of the stationary observer. Similarly, the quantum inequality for the 
flux traveling in one direction is 

We emphasize that these results are only for static two-dimensional spacetimes. The original 
derivation of the quantum inequality of the preceding chapter relies on the metric having a natural 
choice of a timelike Killing vector field. This allows us to clearly define the notion of positive 
frequency for the mode functions. The conformal property of two-dimensional spacetimes is more 
general than the quantum inequality, in that even spacetimes which are time evolving will be 
conformal to some part of two-dimensional Minkowski spacetime. 



Recently, Flanagan [ 64 1 has demonstrated another interesting aspect of the quantum inequalities 
in two-dimensional spacetimes. He has shown that a general form of the quantum inequality can 
be developed in two-dimensional Minkowski spacetime which is independent of the choice of the 
sampling function. We will summarize his results below, but first we will look at the conformal 
triviality of all static two-dimensional spacetimes. 

4.1.1 Conformal Properties 

It is a unique property of two-dimensional spacetimes that any metric can be cast into the form 

g,iu{x) Qfiuiy) = ^'^{y)Vfj,iy (4-3) 

by choosing an appropriate coordinate transform, x ^ y = y{x). For the massless scalar field, the 
wave equation 

U(t){x) = (4.4) 
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then becomes 



n-\y)rj^'''d^d,^{y) = 0. (4.5) 



The mode functions in the two different coordinates are related by 

c^{x) = Ky{x)). (4.6) 

However, (t){y) is the scalar field solution to the wave equation in Minkowski spacetime, composed 
from the standard mode functions 

/fc(y) = (47ru;)-V2e^fc-y, (4.7) 

where ko = oj = \ki\. In two dimensions, the Green's functions are found to transform as 

D{x,x')=D{y{x),y{x')), (4.8) 

where D{y,y') is the standard Minkowski Green's functions. For example, the Wightman Green's 
function in two-dimensional Minkowski spacetime is 

D+{y,y') = {0\^{y)^iy'm = - / -e^Hv-y)^ (4.9) 

47r J -00 uj 

This would also be the Wightman Green's function in an arbitrary curved two-dimensional space- 
time expressed in the coordinates chosen from the metric (|4.3| ). It is evident from the expression 
above that the massless Green's functions have an infrared divergence. It is possible to carry out the 
integration of the above Green's function, and renormalize away the infrared divergence. However, 
for our purposes, the integral representation is sufficient. 

We have already shown for static spacetimes that the QI is given by 

Ap > -^DEDEix, -ito;x, +ito). (4.10) 

If we now make use of the conformal relations for the box operator and the Green's function, we 
find 

Ap> --n-\y)DEDE{y,-ito;y,+ito) = -—^-\y)dl / -e"^-*". (4.11) 

We see that while the Green's function has an infrared divergence, the quantum inequality does 
not. We can interchange the order of the integration and differentiation with the resulting integral 
being well defined for all values of k, and find 

1 

However, in static spacetimes the coordinate time and the proper time of an observer are related 
by To = ^{y)tQ. We can then write the quantum inequality in a more covariant form as 

This is exactly the form of the quantum inequality found by Ford and Roman in Minkowski space- 
time The difference inequality above is applicable to all static two-dimensional spacetimes. 
However the renormalized quantum inequalities in any two-dimensional spacetimes will be different 
because the vacuum energies in these spacetimes are not identical. We can write the renormalized 
two-dimensional quantum inequality as 

PRen. ^ 9 ~l~ Pvacuumi ^^-^^ 

where Pvacuum is the vacuum energy that is defined with respect to the timelike Killing vector field. 
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4.1.2 Sampling Functions in Two-Dimensional Spacetimes 

In the previous chapter, we have proven quantum inequahties using the Lorentzian samphng func- 
tion, Eq. ( |3.12| ). The Lorentzian was chosen to facihtate the proof of the inequahty. Recently, 
Flanagan [^] has shown that quantum inequalities for the massless scalar field in two-dimensional 
Minkowski spacetime can be developed for arbitrary sampling functions. (Note: Flanagan uses the 
term "smearing functions.") We summarize his results below. 

We begin with a smooth, nonnegative function f{v) with unit normalization, 

f{v)dv = l. (4.15) 

It is not necessary that the sampling function have compact support so long as the function ap- 
proaches zero sufficiently quickly to ensure convergence of the above integral. Using the stress-tensor 
in the standard (x, t) coordinates, we define three smeared integral quantities: the time-smeared 
energy density 

/oo 
dtf{t){:Tu{0,t) :), (4.16) 
-oo 

the time-smeared flux 

/oo 
dtf{t){:T-\0,t) :), (4.17) 
-oo 

and the spatially-smeared energy density 

/oo 
dxf{x){:Tu{x,0):). (4.18) 
-oo 

Flanagan then shows that all three of the above integral quantities are bounded below by 

£T,min = £s,min = '^£F,min = —7rr~ / . (4.19) 



247r y_oo f{v) 

This is achieved by use of a Bogolubov transformation which converts the quadratic form of the 



integral in Eq. (4.16|) into a simple form. (Details of the proof can be found in pM-) Flanagan 



also shows that this is an optimum lower bound for a given sampling function. This lower bound is 
achieved by the vacuum state associated with the new coordinates used in the Bogolubov transform. 
This new vacuum state is a generalized multi-mode squeezed state of the original vacuum defined 
in the (x, t) coordinates. 

For arbitrary sampling functions, it is now straightforward to find the optimum quantum in- 
equality in two-dimensional Minkowski spacetime. For example, if we use the Lorentzian sampling 
function, Eq. ( |3.12| ), in the time-smeared energy density, we find 



1 



Ap>-- ^. (4.20) 



This is six times more restrictive than the quantum inequality derived in the preceding section. 
We knew that the quantum inequality developed using the Green's function was a lower bound, 
but it was not necessarily known if it was the optimum lower bound. We now see that Eq. ( 14201 ) 
represents the optimum quantum inequality for the Lorentzian sampling function. 
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4.2 Three-Dimensional Spacetimes 

We now discuss examples of three-dimensional spacetimes in which it is straightforward to find the 
quantum inequalities. First we will develop the quantum inequality in Minkowski spacetime for a 
massive scalar field. The addition of the mass makes the quantum inequality more restrictive. In 
order to generate negative energy densities, we must overcome both the momentum of the particle 
and its rest mass. We will then look at the three-dimensional equivalent of the static Einstein 
universe. Unlike the two-dimensional examples of the preceding section, we will see that in three 
dimensions the form of the quantum inequality is modified from its flat space form because of the 
spacetime curvature. In addition, we will be able to renormalize the quantum inequality in the 
three-dimensional closed universe. Here there exists a non-zero vacuum energy that must be taken 
into consideration. In this case the quantum inequality tells us how much negative energy we can 
measure relative to the background vacuum energy. 

4.2.1 Minkowski Spacetime 

We begin with the easiest example, three-dimensional flat spacetime, where we know that the 
positive frequency massive mode functions are given by 

/k(x, y, t) = [2u;(2^)2]~i/2e.k.x-.c.t ^ (4_21) 



with a frequency to = y/\k\^ + fj? and —oo<ki<oo. The required two-point function is then 
defined by 

1 



/oo 
dk:j: dXy 
-oo 



2 



(4.22) 



v/2a;(27r)2 

It is rather straightforward to carry out the required integrations to find 

GE(2to) = ^e-2-*o. (4.23) 
ovrio 

Now the difference inequality is given by Eq. ( p.23| ), where the Euclidean box operator is 

aE = dl+dl + d^y. (4.24) 

One thing that should be pointed out before proceeding is that the vacuum energy in three- 
dimensional Minkowski spacetime is defined to be zero. Under such circumstances, we are placing 
a lower bound on the energy density of the state itself, and not just on the difference of the energy 
density between the state and the vacuum energy. The quantum inequality becomes 



to 

which upon taking the time derivatives can be written as 



(4.25) 



/5>--^5(2^to). (4.26) 



The "scale function," S{x), is defined as 



S{x) = ( 1 + X + ) e"^ . (4.27) 
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Figure 4.1: The scale function S{2^tQ) for three-dimensional Minkowski spacetime. 



The scale function represents modifications to the quantum inequality due to the mass of the the 
particle. It is plotted in Figure 4.1. When ^ = 0, the scale function becomes one and the quantum 



inequality reduces to the massless form 



P> 



1 



(4.28) 



This is what would have been found if we had chosen the massless mode functions at the beginning 
of this derivation. It is interesting to note that in the limit of the mass of the scalar particle 
becoming very large, the scale function will tend to zero. Thus, if the particle has a mass, it 
becomes more difficult to observe negative energy densities. 

Since the mode functions are invariant under Lorentz transformations, we can always choose 
to develop the inequality in the observer's rest frame and then transform to any other frame. The 
quantum inequality can therefore be written in a more covariant form in flat spacetime as 



■To 

p = — 
IT 



-dr > 



167rT^ 



5(2/iro). 



Here r is the observer's proper time, tq is the sampling time, and u'^ = dx'^ 
four-velocity. 



(4.29) 

{T)/dT is the observer's 



4.2.2 The Closed Universe 

Let us consider the three-dimensional spacetime with a length element given by 

ds^ = -dt^ + (dO'^ + sin^ Od^p^') . 



(4.30) 
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Here constant time-slices of this universe are two spheres of radius a. The wave equation on this 
background with a couphng of strength ^ to the Ricci scalar R is 



- (/u^ + ^R)(f) = 0. 
For the metric ( [4.30 ), the wave equation becomes 



(4.31) 



sm 



{sin 9de<t>) + ^ ^ 2 A "^ + C f^V = 0, (4.32) 
im8 a^'sm^fc' ^ V« / 



which has the solutions 



flmit,0,(p) 



1 



=Y;^(e,(^)e-*"'*. 



(4.33) 



The Yim are the usual spherical harmonics with unit normalization. The mode labels take the 
values, / = 0, 1, 2, • • • and —l<m<l. For each state of the primary quantum number /, there is a 
degeneracy of 2/ + 1 associated with the m, each having the same eigenfrequency 



uJi = a 



V^(^ + 1) + 2e + ^2^2 



(4.34) 



The coupling parameter ^ here can be seen to contribute to the wave functions as a term of the 
same form as the mass. However, we will look only at minimal coupling (,^ = 0). The lower bound 
on the energy density is given by 



1 



00 +1 



2~2uJito 



1=0 m=-l 

However the spherical harmonics obey a sum rule 

+1 



1 CO +1 

^'V^E E -\^Ue,v)\^e~^-^'\ (4.35) 

1=0 m=-l 



8a2 



E V^UOMt 



m=—l 



21 + 1 



(4.36) 



We immediately see that the difference inequality is independent of position, as expected from the 
spatial homogeneity, and that the second term of Eq. (4.35) does not contribute. We have 



A/5 > 



1 



87ra2 



Y^{2l + l)uJie 



2a;; to 



(4.37) 



1=0 



This summation is finite due to the exponentially decaying term. We are now left with evaluating 
the sum for a particular set of values for the mass ^, the radius a, and the sampling time to- Let 
■q = to/a. Then 



Ap > 



1 



2??3^(2Z + l)cD/ e-^^"^' 
1=1 



Fir],p), 



where 



/(/ + 1) +aV- 



(4.38) 



(4.39) 



The coefficient — l/(167rtQ) is the right-hand side of the inequality for the case of a massless field in 
an infinite three-dimensional Minkowski space. The "scale function," F{r],p,), represents how the 
mass and the curvature of the closed spacetime affects the difference inequality. 
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4.2.3 Massless Case 

In terms of the variable rj = to/ a, which is the ratio of the sampling time to the radius of the 
universe, we can write the above expression for F{i], 0) as 



F{r]) = F{i], 0) = 2r/3 ^.{21 + + 1) e^^^^^. (4.40) 

1=1 

A plot of F(r]) is shown in Figure |4.2| . In the limit of ^ 0, when the sampling time is very small 
or the radius of the universe has become so large that it approximates flat space, the function F{rj) 
approaches one, yielding the flat space inequality. 

We can look at the inequality in the two asymptotic regimes of r/. In the large r] regime each 
term of greater / in the exponent will decay faster than the previous term. The / = 1 term yields a 
good approximation. In the other regime, when the sampling time is small compared to the radius, 
we can use the Plana summation formula to calculate the summation explicitly: 



E f{n) + \fm = f{x)dx + r ^^^^^^Sr^dx , (4.41) 



where 

fix) = {2x + l)y'x(x + l)e-W^(^+i). (4.42) 

Immediately we see that for our summation /(O) = 0. The first integral can be done with relative 
ease yielding 

r fix)dx= r(2x + l)Jx(x + 1) e-2W^(^+i)dx = . (4.43) 
Jo Jo ^ 2rj-^ 

This term reproduces the flat space inequality. The second integral in Eq. ( [1.41 ) therefore contains 
all the corrections due to non-zero curvature of the spacetime. Since r] is small, we can expand the 
exponent in a Taylor series around r] = 0. Keeping the lowest order terms, we find 

^^'''^ - ^^-'^Idx ~ -/o - r?/i + • • • , (4.44) 



o2lTX ]^ 



where 



and 



V2^ ^ ^ dx PS 0.265096, (4.45) 



2x\l ^Jx"^ + 1 — X + y \/x^~+T + x 
e^'^^ - 1 

4:x{2x^ - 1) , 2 



dx = -—. (4.46) 



Therefore the function F{r]) in the small rj limit is given by 

F(r]) ~ 1 - 0.530192 rf + — + 0(rf) + • • • , (4.47) 

15 

and in the large r] limit by 

F(r/) ~ QV2rf e'^^'' + • • • . (4.48) 



Both of these asymptotic forms are plotted along with the exact form of F{ri) in Figure [4.2[ The 
graph shows that the two asymptotic limits follow the exact graph very precisely except in the 
interval 1 < r/ < 2. These results for the function F{r]), combined with Eq. (f4.38| ), yield the 
difference inequality for the massless scalar field. 



CHAPTER 4. SCALAR FIELD EXAMPLES 



46 



F(T1) 




1 2 3 4 5 



Figure 4.2: The scale function F{r]) and its asymptotic forms for the three-dimensional closed uni- 
verse. The solid line is the exact form of the function. The dotted line is the small rj approximation, 
while the dot-dash line is the large rj approximation. 
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The difference inequality does not require knowledge of the actual value of the renormalized 
vacuum energy, Pvacuum- However, if we wish to obtain a bound on the energy density itself, 
we must combine the difference inequality and Pvacuum- One procedure for computing Pvacuum is 
analogous to that used to find the Casimir energy in flat spacetime with boundaries: One defines 
a regularized energy density, subtracts the corresponding flat space energy density, and then takes 
the limit in which the regulator is removed. A possible choice of regulator is to insert a cutoff 
function, g{u))^ in the mode sum and define the regularized energy density as 

1 -j^ oo 

Preg = ^ E sI^a) = ^ + 1) UJi g{uJi) . (4.49) 

A 1=0 

In the limit that a ^ oo, we may replace the sum by an integral and obtain the regularized flat 
space energy density: 

1 

PFSreg = -r dujujg{uj). (4.50) 
47r JO 

The renormalized vacuum energy density may then be defined as 

Pvacuum = ^^^{Preg — PFSreg) ■ (4-51) 

The vacuum energy density so obtained will be denoted by pcasimir- An analogous procedure 
was used in |jl^ to obtain the vacuum energy density for the conformal scalar field in the four- 
dimensional Einstein universe. An explicit calculation for the present case, again using the Plana 
summation formula, yields 

- _ ~ 0-265096 

PCasimir — o ^ o ■ V^-^^) 



This agrees with the result obtained by Elizalde [15| using the zeta function technique. In this 
case PCasimir < 0, whereas the analogous calculation for the conformal scalar field in this three- 
dimensional spacetime yields a positive vacuum energy density, pcasimir — 

l/967ra3. It should 

be noted that this procedure works for the closed universe because the divergent part of preg is 
independent of a. More generally, there may be curvature-dependent divergences which must also 
be removed. 

Let us now return to the explicit forms of the difference inequality. In the limit that Iq ^ a, 
Eqs. (l438|) and (gl|) yield 

Ap>-^e-2^WV (453) 
Similarly, in the limit that we find 

^ - 16^ - - 6^ + • • • ■ ^^-^^^ 

Thus the bound on the renormalized energy density in an arbitrary quantum state in the latter 
limit becomes 

In the a ^ oo limit, both Eqs. ([4.54| ) and ( [4.55| ) reduce to the quantum inequality for three- 



dimensional Minkowski space, Eq. (|]2|). In this limit 

PCasimir — SO — p. In the other 



limit when the sampling time becomes long, we find that Ap decays exponentially as a function 
of the sampling time. This simply reflects the fact that the difference in energy density between 
an arbitrary state and the vacuum state satisfies the quantum averaged weak energy condition, 
Eq. (|L1|). 
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4.3 Four-Dimensional Robertson- Walker Universe 

Now we will apply the same method to the case of the three homogeneous and isotropic universes 
given by the four-dimensional static Robertson- Walker metrics. Here we hav^ 

[e = 0] : gijdx'dx^ = dx^ + dy'^ + dz^, (4.56) 

for the flat universe with no curvature (Minkowski spacetime). For the closed universe with constant 
radius a, i.e., the universe of constant positive curvature, the spatial length element is given by 

[e = 1] : g^jdx'dx^ = \dx^ + sin^ xidO"^ + sin^ Od^"^)] , (4.57) 

where 0<x^7r, O<0<7r, and < 99 < 2tt. The open universe [e = —1] is given by making the 
replacement sinx sinhx in Eq. ( 4.57| ) and now allowing x to take on the values < x < c«. To 



find the lower bound of the quantum inequalities above we must solve for the eigenfunctions of the 
covariant Helmholtz equation 

ViV^C/A(x) + {ujI - fi^)Ux{^) = 0, (4.58) 

where // is the mass of the scalar field and uj\ is the energy. A useful form of the solutions for this 
case is given by Parker and Fulling [p^. 

4.3.1 Flat and Open Universes in Four Dimensions 

In the notation developed in Chapter |3|, the spatial portion of the wave functions for flat (Euclidean) 
space is given by (continuum normalization) 

[e = 0] : ?7k(x) = [2u;(27r)3]-V2gik.x^ (4_5g) 



to = ^|k|2 + /i2^ (460) 

k = (fci, A;2, /C3) (— 00 < A;j < 00). 

It is evident that |C/kP will be independent of position. This immediately removes the second term 
of the inequality in Eq. ( p.l8[ ). 



In the open universe, the spatial functions are given by 

[e = -l]: C/a(x) = {2a'u;,)-'/'u';;;\x)Yim{0,'p), (4.61) 

A = {q,l,m). 

Here < g < 00; / = 0, 1, • • • ; and m = —I, —l + l,---, +1. The sum over all modes involves an 
integral over the radial momentum q. The functions H^^ \x) ci-re given in Eq. (5.23) of Birrell and 
Davies |53|. Apart from the normalization factor, they are 



' (X) oc sinh' X (^ ^cQgi^^ j cos qx- (4.63) 



^We are using e as a free parameter which takes the value e = — 1 in the open universe, e = in the flat universe, 
and e = +1 in the closed universes. 
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As with the mode functions of the three-dimensional closed spacetime above, the mode functions 
of the open four-dimensional universe satisfy an addition theorem [ 



67|, |68| 



Im 



27r2 • 



(4.64) 



Since the addition theorem removes any spatial dependence, we again get no contribution from the 
Laplacian term of the quantum inequality, Eq. ( p. 18 ). Upon substitution of the mode functions for 
both the flat and open universes into the quantum inequality and using the addition theorem in 
the open spacetime case we have 



and 



0] 



-1] 



Ap > 



1 



167r3 7_ 



d^kuj]f^ e' 



-2a;kto 



1 1'°° 



■2uj„to 



(4.65) 



(4.66) 



respectively. The three-dimensional integral in momentum space can be carried out by making a 
change to spherical momentum coordinates. The two cases can be written compactly as 



dk A:2we-2'^*o, 



where 



u) = Lj{k, a, /i) = (jt"^ — e/a^ + fi^ 



1/2 



(4.67) 



(4.68) 



Note that e = for flat space and k = q/a for the open universe. This integral can be carried out 
explicitly in terms of modified Bessel functions Kn{z). The result is 



Ap > 



327r2t4 



Uz^K,[z)-z''K2{z)) 



G(z), 



where 



2toW/x 



i2 • 



(4.69) 



(4.70) 



The coefficient — 3/(327r2tQ) is the lower bound on p found in 48| for a massless scalar field 
in Minkowski spacetime. The function G{z) is the "scale function," similar to that found in the 
case of the three-dimensional closed universe. It is the same function found in for Minkowski 
spacetime (e = 0), and is plotted in Figure 4^. Again we see that in the limit of z ^ the 
scale function approaches unity, returning the flat space massless inequality in four dimensions, 
Eq. (ra . 



4.3.2 The Einstein Universe 

In the case of the closed universe, the spatial mode functions are the four-dimensional spherical 
harmonics, which have the form |6(;, ^ 

[€ = !]: C/a(x) 



(2ai,a3)-i/2nW(;t)Y;™(^,<^), 



n{n + 2) 
(n, /, m). 



(4.71) 
(4.72) 
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Here, n = 0, 1, 2, • • • ; / = 0, 1, • • • , n; and m = —I, — / + 1, • • • , +1. The function II^"'") is found from 
n^"^ by replacing x by —ix and q by —i{n + 1) |66]. Alternatively, they can be written in terms of 
Gegenbauer polynomials 111, ^ as 



n(;[)(x)ocsin'xC^+_'/(cosx). 



(4.73) 



In either case, the addition theorem is found from Eq. 11.4(3) of [|69| ] (with p = 2 and = r]). This 
reduces to 



J2\KVixmmie,^)\' 



Im 



2^2 



from which it is easy to show that the energy density inequality, Eq. ( p.l8| ), becomes 

1 



47r2a3 



n=0 



If we use the variable rj = to/a in the above equation, we can simplify it to 

3 



A/5 > 



327r2f4 



(4.74) 



(4.75) 



(4.76) 



Again we find the flat space solution in four dimensions, multiplied by the scale function H{r},fj,), 
which is defined as 



H{ri, ^) = -7]^ ^ (n + 1) Vn(n + 2) + a'^fi'^ e 



,2, ,2 -2?7^n(n+2)+a2^2 



(4.77) 



n=l 



and is plotted in Figure i.4 for fi = 0. The scale function here has a small bump occurring at 



roughly rj ~ 0.5 with a height of 1.03245. This may permit the magnitude of the negative energy 
to be slightly greater for a massless scalar field in the Einstein universe than is allowed in a flat 
universe for comparable sampling times. A similar result was shown to exist for massive fields in 



two-dimensional Minkowski spacetime |48]. 



4.3.3 Massless Asymptotic Limits in the Einstein Universe 

As with the three-dimensional closed universe, we can find the asymptotic limits of this function. 
We again follow the method of the previous section, assuming the scalar field is massless, and 
making use of the Plana summation formula to find 



where 



and 



I2 = J^{x + l)yx(x + 2) e" W^(^+2) dx , 



2Re 



(ix + ifJixiix + 2)e-2''V«^(«+2) 

I I — dx 





327rx 



2x - 


-0.356109. 



l)\J^/x^ + A + X- 2x\Jy/x'^ + 4 
e2'^^ - 1 



dx , 



(4.78) 
(4.79) 



(4.80) 
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The first integral can be done in terms of Struve H„(2) and Neumann Nn{z) functions, with the 
result 



16 dr}^ 



;Hi(2r?)-iVi(2r/)) 



4 

H f] ■ 

15 ' 



(4.81) 



If we follow the same procedure as in the previous section for defining the renormalized vacuum 
energy density for the minimally coupled scalar field in the Einstein universe, then we obtain 



1 



PCa 



0.356109 



a" 



l/4807r^a^ for the massless conformal scalar field |13 
differs from that obtained by Elizalde 



(4.82) 



Here 



The same method yields pca 

our result for the minimal field, Eq. ( 4.82 ), differs from that obtained by Elizalde |[5[| using the 
zeta function method, p'casimir ~ — 0.411502/47r^a^. This discrepancy probably reflects that the 
renormalized vacuum energy density is not uniquely defined. The renormalized stress-tensor in 
a curved spacetime is only defined up to additional finite renormalizations of the form of those 
required to remove the infinities. In general this includes the geometrical tensors ^^^H^y and 



55, 70 1 for the definitions of these tensors and a discussion 



L^y. (See any of the references in [ 
of their role in renormalization.) In the Einstein universe, both of these tensors are nonzero and 
are proportional to 1/a^. Thus the addition of these tensors to {T^^) will change the numerical 
coefficient in pcasimir- The logarithmically divergent parts of {T^y) which are proportional to 
^^^H^y and ^"^^H^^ happen to vanish in the Einstein universe, but not in a more general spacetime. 
In principle, we should imagine that the renormalization procedure is performed in an arbitrary 
spacetime, and only later do we specialize to a specific metric. Unfortunately, it is computationally 
impossible to do this explicitly. Thus, the fact that a particular divergent term happens to vanish 
in a particular spacetime does not preclude the presence of finite terms of the same form. 

In the small r] limit, we can expand the Struve and the Neumann functions in Eq. ( [4.81| ) in a 
Taylor series to obtain 



1 



//(??, 0) = l + -r/^ + 



(3 + 47) , 8 



1 



r + -(-0.356109)??* + -r?* Inr? + 0{rf) + 



(4.83) 



3' 12 ' 3' "3 

where 7 is Euler's constant, which arises in the Taylor series expansion of the Neumann function. 
This is similar to that of the three-dimensional universe and again contains a term of the form of 



the Casimir energy. In the large rj limit we again keep just the first term of the series (4.77). Both 



asymptotic forms are plotted with the exact solution in Figure [4.4| . We see that the asymptotic 
form is again a very good approximation except in the interval 1 < r/ < 2, as was the case for the 
three-dimensional closed universe. 

The difference inequalities for massless fields are then given by 



327r2t4 



1 2 , (3 + 47) 4^14, ^ 
3^ 12 — ^ 3^ 



PCasimir 



for Iq <^ a and 



%/3 



for a . 



(4.84) 



(4.85) 



Using A/5 = pRen. - PCasimir-, where PCasimir is the expectation value of pRen. in the vacuum 
state, we can calculate the renormalized energy density that would be constrained by the quantum 
inequalities, subject to renormalization ambiguities. For example, in the four-dimensional Einstein 
universe the renormalized energy density inequality is 



PRe 



> 



32^2^4 



^ 1 2 , (3 + 47) 4 1 4 

1 + -ri + ^ -'n + -fj 

3 ' 12 ' 3 ' 



Inr? + 0{r]'^) + • 



(4.86) 
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Figure 4.4: The scale function H{r]) for the Einstein universe. The sohd hne is the exact result. 
The dotted line is the asymptotic expansion for small ij, while the dot-dash curve is the large r] 
approximation. The maximum occurs at « 0.5. 

for Iq <^ a. Here the coefficient of the r/^ term could be altered by a finite renormalization, but 
the rest of the expression is unambiguous. This result is identical to the short sampling time 
expansion of the quantum inequality in the Einstein universe, Eq. ( 3.57| ), developed in Section |3.3| . 
Expressions similar to Eq. ([4. 861) could be found for any of the other cases above. In the case of 



the flat Robertson- Walker universe, there is no Casimir vacuum energy. Under such circumstances 
the difference inequality and the renormalized energy density inequality are the same, and are free 
of ambiguities. 



4.4 Quantum Inequalities Near Planar Mirrors 
4.4.1 Single Mirror 

Consider four-dimensional Minkowski spacetime with a perfectly reflecting boundary at z = 0, 
located in the x — y plane, at which we require the scalar field to vanish. The two-point function 
can be found by using the standard Feynman Green's function in Minkowski space. 



{X,X )- _ ^,)2 + (y _ y/)2 + _ ^/)2 _ _ i/)2] ' 



(4.87) 



and applying the method of images to find the required Green's function when the boundary is 
present. For a single conducting plate we have 

1 



Gf{x,x') = ^ 



{x - x'Y + (y - y'Y + (z - z'Y - (i - t'f 
1 



{x - x'Y + (y - y'Y + (z + z'Y -{t- t'Y 
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If we Euclideanize by allowing t — > —Uq, t' — > Uq and then take x' x, we find 

GE{2to) = T^(i- ,A-o) . (4.89) 



In addition, the Euclidean box operator is given by 

aE = dl + dl + d^ + dl (4.90) 

It is easily shown that the quantum inequality is given by 

Ap>4n.G.(2t„) = -^ + :j^j-j^i^. (4.91) 

The first term of this inequality is identical to that for Minkowski space. The second term 
represents the effect of the mirror on the quantum inequality. For the minimally coupled scalar field 
we know from Eq. ( |2.66| ) that there is a non-zero, negative vacuum energy density which diverges 
as the mirror is approached. Adding this vacuum term to both the left- and right-hand sides of the 
above expression allows us to find the renormalized quantum inequality for this spacetime, 

3 11 , , 



32tt^4 167r2(t2 + ^2)2 167^2 



There are two limits in which the behavior of the renormalized quantum inequality can be 
studied. First consider z ^ tQ. In this limit, the correction term due to the mirror, and the 
vacuum energy very nearly cancel and the quantum inequality reduces to 

This is exactly the expression for the quantum inequality in Minkowski spacetime. Thus, if an 
observer samples the energy density on time scales which are small compared to the light travel 
time to the boundary, then the Minkowski space quantum inequality is a good approximation. 

The other important limit is when z <^ Iq. This is the case for observations made very close to 
the mirror, but for very long times. The quantum inequality then reduces to 

PKen.>-^. (4.94) 

Here, we see that the quantum field is satisfying the quantum averaged weak energy condition. 
Recall that throughout the present paper, we are concerned with observers at rest with respect to 
the plate. If the observer is moving and passes through the plate, then it is necessary to reformulate 
the quantum inequalities in terms of sampling functions with compact support [^l|. It should be 
noted that the divergence of the vacuum energy on the plate is due to the unphysical nature of 
perfectly reflecting boundary conditions. If the mirror becomes transparent at high frequencies, 
the divergence is removed. Even if the mirror is perfectly reflecting, but has a nonzero position 



uncertainty, the divergence is also removed |72]. 



4.4.2 Two Parallel Plates 

Now let us consider the case of two parallel plates, one located in the z = plane and another 
located in the z = L plane. We are interested in finding the quantum inequality in the region 
between the two plates, namely < z < L. We can again use the method of images to find the 



CHAPTER 4. SCALAR FIELD EXAMPLES 



55 



Green's function. In this case, not only do we have to consider the reflection of the source in 
each mirror, but we must also take into account the reflection of one image in the other mirror, 
and then the reflection of the reflections. This leads to an infinite number of terms that must be 
summed to find the exact form of the Green's function. If we place a source at {t' ,x' ,y' , z'), then 
there is an image of the source at {t',x',y', —z') from the mirror at z = and a second image at 
{t' , x', y', 2L — z') from the mirror at z = L. Then, we must add the images of these images to the 
Green's function, continuing ad infinitum for every pair of resulting images. If we use the notation 

Gf{z, a± z') = Gp^ (t, X, y, z; t' ,x' , y' , a ± z'), (4.95) 

where a is a constant, then we can write the Green's functions between the plates as 

oo 

Gix,x') = GFiz,z')-GFiz,-z') + J2[GFiz,2nL + z')-GFiz,-2nL-z') 

n=l 

+Gf(z, -2nL + z') - Gf{z, 2nL - z')] . (4.96) 

We Euclideanize as above, and let the spatial separation between the source and observer points 
go to zero, to find 



le^Hig + 



oo , 2 1 1 



167r2 ^ 

n=l 



tl + (nL)2 tl + (nL + zf 1% + [nl 
It is now straightforward to find the quantum inequality, 



(4.97) 



3 1 

Ap > -— ^ + 



327r2t4 167r2(i2 + ^2)2 

1 ^ UuLf-^tl X X 1 

^ 1 r+2 I r ^213 ^ [+2 I r, _L ^^212 ^ [+2 , r _ ^^212 ( ' V*-"^"; 



167r2 \ {tl + (nL)2]3 [tg + (^^ + ^)2]2 [^2 + (^^ _ ^)2p 

The first term in the expression above is identical to that found for Minkowski space. The second 
term is the modification of the quantum inequality due to the mirror at z = 0. The modification 
due to the presence of the second mirror is contained in the summation, as well as all of the multiple 



reflection contributions. When the Casimir vacuum energy, given by |55] 



__^ 3-2sin2(../L) 
Pvacuum- ^^^4 sin4(7rz/L) 1440L4 ' ^^-^^^ 

is added back into this equation for renormalization, we find, as with a single mirror, that close 
to the mirror surfaces the vacuum energy comes to dominate and the quantum inequality becomes 
extremely weak. 

4.5 Spacetimes with Horizons 

We will now change from fiat spacetimes with boundaries to spacetimes in which there exist hori- 
zons. We will begin with the two-dimensional Rindler spacetime to develop the quantum inequality 
for uniformly accelerating observers. For these observers, there exists a particle horizon along the 



null rays x = =bt (see Figure 4.5). 
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We will then look at the static coordinatization of de Sitter spacetime. Again there exists a 
particle horizon in this spacetime, somewhat similar to that of the Rindler spacetime. The two 
problems differ because Rindler spacetime is flat while the de Sitter spacetime has constant, positive 
spacetime curvature. 



4.5.1 Two-Dimensional Rindler Spacetime 

We begin with the usual two-dimensional Minkowski metric 

ds^ = -dt^ + dx^ . (4.100) 

Now let us consider an observer who is moving with constant acceleration. We can transform to 
the observer's rest frame (Section 4.5 of [^]) by 

t = a'^e^^sinhar?, (4.101) 
X = o~^e"^cosha?7, (4.102) 

where a is a constant related to the acceleration by 

a e~°'^ = proper acceleration. (4.103) 

The metric in the rest frame of the observer is then given by 

ds"^ = e'^<{-dif + di^). (4.104) 

The accelerating observer's coordinates (?7, ^) only cover one quadrant of Minkowski spacetime, 
where x > \t\. This is shown in Figure |4.5| . Four different coordinate patches are required to cover 
all of Minkowski spacetime in the regions labeled L, R, F, and P. For the remainder of the paper we 
will be working specifically in the left and right regions, labeled L and R respectively. In these two 
regions, uniformly accelerating observers in Minkowski spacetime can be represented by observers 



at rest at constant ^ in Rindler coordinates, as shown by the hyperbola in Figure |4.5 . 
The massless scalar wave equation in Rindler spacetime is given by 

which has the positive frequency mode function solutions 

A(r/,0 = (47ru;)-^/2e**=^±*-''. (4.106) 

Here — oo < k < oo and to = \k\. The plus and minus signs correspond to the left and right Rindler 
wedges, respectively. Using the above mode functions, we can expand the general solution as 

HV,0 = r dk \biMv,0+bj:^rk{v,0+b§h{v,0+b§^fkiv,0] , (4.107) 

where and 6^ are the creation and annihilation operators in the left Rindler wedge and similarly 

for b^^ and bj} in the right Rindler wedge. We also need to define two vacua, |0l) and |0_r), with 
the properties 

bk^m = b§^\OL) = b^\OL) = bj^\OR) = b§\OL) = bj:\OR) = 0. (4.108) 
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Figure 4.5: The Rindler coordinatization of two-dimensional Minkowski spacetime. The time co- 
ordinate rj = constant are straight Hnes passing through the origin, while the space coordinate ^ = 
constant are hyperbolae. The Minkowski spacetime is covered by four separate coordinate patches, 
labeled by L, R, F, and P. The two null rays {x = t and x = —t) act as horizons. 
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The Rindler particle states are then excitations above the vacuum given by 

|{U}l> = bi^\OL), (4.109) 
= b§^\On). (4.110) 

With this in hand, we can find the two-point function in either the left or right regions. Let us 
consider the right Rindler wedge, where 

G+ix,x') = {OR\<P{x)<Pix')\OR), 
dkfk{x)fl{x) , 



= J_ f°° ^gifc(€-€')-*^{'?-'?') (4 111) 

To find the Euclidean two-point function required for the quantum inequality, we first allow the 
spatial separation to go to zero and then take {rj — rj') — 2ir/o, yielding 

G^(2%) = ^ (4.112) 
2tt Jo lu 

In two dimensions, the Euclidean Green's function for the massless scalar field has an infrared 
divergence as can be seen from the form above, in which the integral is not well defined in the 
limit of u; ^ 0. However, in the process of finding the quantum inequality we act on the Green's 
function with the Euclidean box operator. If we first take the derivatives of the Green's function, 
and then carry out the integration, the result is well defined for all values of uj. In Rindler space, 
the Euclidean box operator is given by 



°E = e-^< + . (4.113) 




It is now easy to solve for the quantum inequality 

1 1 r°° 1 

> —UEGE{2m) = -T^e-2< / da; u;e~^^^'' = ——^ . (4.114) 

4 ^ ' 2t: Jo 8^(e<r/o) 

However, the coordinate time r]Q is related to the observer's proper time by 

TO = e<r/o , (4.115) 

allowing us to rewrite the quantum inequality in a more covariant form, 

1 

^0 



Ap> (4.116) 



This is exactly the same form of the quantum inequality as found in two-dimensional Minkowski 
spacetime [32, 48|. As we have seen in Section [4.1.1| , this is a typical property of static two- 



dimensional spacetimes which arises because all two-dimensional static spacetimes are conformal 
to one another. However, we reiterate that the renormalized quantum inequalities are not identical 
in different spacetimes because of differences in the vacuum energies. 
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4.5.2 de Sitter Spacetime 

Let us now consider four-dimensional de Sitter spacetime. The scalar field quantum inequality, 
Eq. ( |3.23| ), assumes a timelike Killing vector, so it will be convenient to use the static parametriza- 
tion of de Sitter space, 



-(1 - ^) dt^ + (1 - ^)~^dr^ + r^{de^ + sin^ 6 dip'') . 



a" 



(4.117) 



There is a particle horizon at r = a for an observer sitting at rest at r = 0. The coordinates take 
the values, 0<r<a, O<0<7r, and < 99 < 27r. It should be noted that this choice of metric 
covers one quarter of de Sitter spacetime. 
The scalar wave equation is 



1-^1 d^-^d.^ 



l-^\dr 



a" 



where is the square of the standard angular momentum operator and is defined by 



1 



1 



-de{siTiede) + —^dl. 
smy sm''^ ^ 



(4.118) 
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The unit norm positive frequency mode functions are found [ |73| , [f^ , |75| , 76, 77] to be of the form 

1 



4>uj,l,ni{t,r,9,(p) 



-iujt 



(4.120) 



UJ 



where z = r/a is a dimensionless length, Yim{0,ip) are the standard spherical harmonics, and the 
mode labels I and m take the values ^ = 0, 1, 2, • • • and —l<m<l. The radial portion of the the 
solution is given by 



T{1 + ^)T{iaoj) 



^ (1 - z2)W2 F(5-^ 5+. I + ^. ^2) 



(4.121) 



where F(q;, /3; 7; 2;) is the hypergeometric function |78] and 



''I 



- I / + - + iauj lb ^ / - — a^/i^ 



(4.122) 



We can then express the two-point function as 

/ dk—^ 

• L 4-71 

Im 



g{x,x')=y: / 



fi*iz) fUz')YUO,^)yim{0',v') e 



1 1 J 



/ ,l\ ik{t-t')/a 



(4.123) 



where k = auj. Now if we Euclideanize according to Eq. ( 3.21 ) and set the spatial separation of the 
points to zero, we may use the addition theorem for the spherical harmonics, Eq. ( [4. 361) , to find 
the Euclidean Green's function 



Ge 



1 



IGvr^a^ 



k 



T{1 + ^)T{tk) 



.21 



(4.124) 



This is independent of the angular coordinates, as expected, because de Sitter spacetime is isotropic. 
We now need the Euclidean box operator. Because of the angular independence of the Green's 
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function, it is only necessary to know the temporal and radial portions of the box operator. The 
energy density inequality, Eq. (l3.23| ), then becomes 



Ap > 



-dl + 



5. 



z^il-z'^)d, \ G'i5(x,-to;x,+to). 



(4.125) 



The temporal derivative term in Eq. ( 4.125 ) will simply bring down two powers of k/a. Using 
the properties of the hypergeometric function, it can be shown that 



(1-zT F\bT,bi;l + l;z% 



(4.126) 



from which we can take the appropriate spatial derivatives. If we allow z ^ 0, then we have F — > 1 
and only the / = terms will contribute in the time derivative part of Eq. (4.125). For the radial 
derivative, it is possible to show that 



hm \z' (1 - z')d, \z'^{l - zy'F\z')] I 



2^U Z 



2(aV^ - k"^) for I = 0, 
6 for / = 1, 

otherwise. 



(4.127) 



Using this result, we find for the observer at r = that 



A/5 > 



1 



dk sinh(7rA;) 



r(6o-)r(6+) +4r(6r)r(6+) 



-2tok/a 



. (4.128) 



There are two cases for which the right-hand side can be evaluated analytically, /i = and /i 
V2/a. For fi = 0, we have 



Ap > 



Svr^a^ Jo 



1 



k. 



Stt'^u'^ Jo 



dk sinh(7r/c) 

) 

dA:(2A:3 + 5A;)e-2*o'=/" 



.k. 



^(^2)^(2+^2) 



+ 4 



1 + 



5 /to 



3 V a 



-2tok/a 



(4.129) 



where we have used the identities for gamma functions proven in Appendix Similarly for p 
\/2/a, we find 



Ap > 



327r2t4 



1 + 



(4.130) 



We can compare these results with the short sampling time approximation from Section Solving 
for the necessary geometrical coefficients, we find 



^000 



vi 



'29 1 l/i^' 
60 c? ~ 4^ 



\9tt\ 



29 J_ _ Ifj^ 1 4 
60^ ~ 2^ ^ 8^ ' 

1 (3-r7a^) 



Q,2 _ ^2 jop.'^ 

The general short sampling time expansion, Eq. ( ^.52] ), now becomes 



(4.131) 

(4.132) 
(4.133) 
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A/5 > 



327r2i 



1 + 



^ 1 (3 - rVa^) 

Q,2 _ ^2 Iq^2^ 



+^ f/^'-A)^oM2roVa^) + OK) + 



(4.134) 



where tq = (l — r"^ /a'^)^^'^to. If r = and /u takes the values or y/2/a, this agrees with Eqs. ( [4.12S| ) 
or ( [4.130| ), respectively. Note that this small to expansion is valid for all radii, < r < a. We can 
also find the proper sampling time from Eq. ( 3.54| ) for which this expansion is valid, 



To <. Tr, 



5 — Sr'^/a'^ 



(4.135) 



For the observer sitting at the origin of the coordinate system, tq <C aj^fh. This is the scale on 
which the spacetime can be considered "locally flat." For observers at r > 0, who do not move on 
geodesies, Tm. decreases and approaches zero as r 



a: 



Tm ~ \la[a - r), 
Note that the proper distance to the horizon from radius r is 
dr' 



a. 



i 



r ^1 - r''^/a^ 



a [7r/2 — arcsin(r/a)] ~ J2a{a 



as r 



a. 



(4.136) 



(4.137) 



Thus, for observers close to the horizon, if the sampling time is small compared to this distance to 
the horizon, i.e., if tq <^ £, then tq <C and the short sampling time expansion is valid. 

We can also obtain a renormalized quantum inequality for the energy density at the origin for 
the case /i = \/2a. By the addition of the vacuum energy to both sides of Eq. ( [4.130| ) we find for 
r = 



> 



327r2t4 



1 + 



a J 



1 



9607r2a4 



(4.138) 



We can now predict what will happen in the infinite sampling time limit of the renormalized 
quantum inequality for any observer's position. We know from Eqs. ( [4.124| ) and ( [4.125| ) that the 
difference inequality will always go to zero, yielding a QAWEC in static de Sitter spacetime of 



lim — 

to^oo vr J_ 



I Ren. 



dt > 



1 



4807r2a^ 



a 



(a2 



1 



(4.139) 



We immediately see that for a static observer who is arbitrarily close to the horizon in de Sitter 
spacetime, the right hand side of Eq. (4.139) becomes extremely negative, and diverges on the 
horizon itself. This is similar to the behavior found for static observers located near the perfectly 
reflecting mirror discussed earlier and for Rindler space. 



4.6 Black Holes 



We now turn our attention to an especially interesting spacetime in which quantum inequalities 
can be developed, the exterior region of a black hole in two and four dimensions. 
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4.6.1 Two-Dimensional Black Holes 

Let us consider the metric 

ds'^ = -C{r)df + C{r)-^dr^ , (4.140) 

where C(r) is a function chosen such that C ^ 1 and dC/dr ^ as r ^ oo. Additionally, there is 
an event horizon at some value tq where C(ro) = 0. For example, in the Schwarzschild spacetime, 
C(r) = 1 - 2Mr"\ there is a horizon at r = 2M. Another choice for C is that of the Reissner- 
Nordstrom black hole, where C(r) = 1 — 2Mr~^ + Q^r~^. In general, we will leave the function C 
unspecified for the remainder of the derivation. The above metric leads to the massless, minimally 
coupled scalar wave equation 

■dl(t){r,t) + dr[C{r)dr<p{r,t)\=Q. (4.141) 



C{r. 



Unlike in four dimensions, the two-dimensional wave equation can be analytically solved everywhere. 
If we use the standard definition of the r* coordinate, 

then it is convenient for us to take as the definition of the positive frequency mode functions 

/fe(r, t) = i (47ra;)-^/2 ^ikr*-i.t^ ^ ^ |^|^ (4 ^43) 

where —00 < A; < 00. 

The problem of finding the quantum inequality simply reduces to using the mode functions to 
find the Euclidean Green's function. We have 

dk ^=e'^^* e-2'"*° = — / dww-V^'^*". (4.144) 
-00 v47ra; 27r 

As in the case of two-dimensional Rindler space, the Euclidean Green's function has an infrared 
divergence. Wc can again apply the Euclidean box operator first and then do the integration to 
obtain the quantum inequality, 

1 r°° 1 

> -— / duj u e-2'^*° = , . . (4.145) 

- 2TrC{r)Jo 8TTC{r)tl ^ ' 

However, the observer's proper time is related to the coordinate time by r = C(r)^/^i, such that 
we can write the difference inequality as 

Ap>-^. (4.146) 

This is the same form as found for two-dimensional Minkowski and Rindler spacetimes. This is the 
expected result because of the conformal equivalence of all two-dimensional spacetimes. 

This now brings us to the matter of renormalization. There exist three candidates for the 
vacuum state of a black hole: the Boulware vacuum, the Hartle-Hawking vacuum, and the Unruh 
vacuum. However the derivation of the difference inequality relies on the mode functions being 
defined to have positive frequency with respect to the timelike Killing vector dt-, and the vacuum 
state being destroyed by the annihilation operator, i.e., 



O'k |Ofc) = 0, for all k . 



(4.147) 
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In Schwarzschild spacetime, this defines the Boulware vacuum. Thus, we can solve for the renor- 
mahzed quantum inequahty, 



PRe 



{Tu/\gtt\) 



Ren. 



+ 



dt > 



1 



87rr2 



+ PB{r) . 



(4.148) 



The Boulware vacuum energy density in two dimensions for the Reissner-Nordstrom black hole is 
given explicitly by (see Section 8.2 of 



247r \ r r'^ 



AM 7M2 6Q2 14MQ2 5g4 



In the limit tq — > oo, we recover a QAWEC condition on the energy density 

{Ttt/gtt) Ren. 



y to 

lim — 



-^T^^dt > PB{r). 



(4.149) 



(4.150) 



This has the interpretation that the integrated energy density in an arbitrary particle state can 
never be more negative than that of the Boulware vacuum state. In particular, this will be true for 
the Hartle-Hawking and Unruh vacuum states. 



4.6.2 Four-Dimensional Schwarzschild Spacetime 

Now let us turn to the four-dimensional Schwarzschild spacetime with the metric 



ds' 



2M^ 

,dr+[l- 
r J \ r J 



2M\ , 9 , 



(4.151) 



The normalized mode functions for a massless scalar field in the exterior region (r > 2M) of 
Schwarzschild spacetime can be written as [54| 



{AiTufl^-'^' Rl (^|r)Y,„(0,(^), 
(47ra;)i/2e-^"* % {Lo\r)Y ^9 , ^) , 



(4.152) 



where Ri {oj\r) and Ri (w|r) are the outgoing and ingoing solutions to the radial portion of the 
wave equation, respectively. Although they cannot be written down analytically, their asymptotic 
forms are 



Rl {^\r) 

for the outgoing modes and 
Rl (wk) 



r e 



+ A; {uj)r 



— 1 „~iu)r* 



+ A; (u;)r-ie''^''* 



2M, 

oo, 

2M, 
oo. 



(4.153) 



(4.154) 



for the ingoing modes. The normalization factors B;(c<j), A; (i^), and Ai (w) are the transmission and 
reflection coefficients for the scalar field with an angular momentum-dependent potential barrier. 
Now let us consider the two-point function in the Boulware vacuum. It is given by 



Gs(x,x') = E 



Im 



doj 

( 

Attuj 



Y,^(0,^)YL(0',^') 



Rl {uj\r) Rl {uj\r')+ Ri {uj\r) Ri {uj\r') . 

(4.155) 
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We are interested in the two-point function when the spatial separation goes to zero, i.e., letting 
r' ^ r, 9' ^ 9, and if' ^ if. We can again make use of the addition theorem for the spherical 
harmonics, Eq. (4.36). Let us also Euclideanize, by taking {t—t') —2itQ. The Euclidean two-point 
function then reduces to 



GBEC^ta 



1 



167r2 



duj 



J2 ^e-2*o(2/ + l) \Ri{u;\r)\^ + \m{u\r) 







(4.156) 



In the two asymptotic regimes, close to the event horizon of the black hole (r 2M), or far from 
the black hole (r — > oo), the radial portion of the wave equation also satisfies a sum rule. It was 
found by Candelas that 



LO r 



and 



> J 4^2(1 - 2M/r)-\ 

(2M)-2E£o(2/ + l)|BKu;)|2, 
4a;2^ 



5:(2/ + l)|ft {co\r)\'r. 
1=0 

with the coefficient B;(a;) given, in the case 2Muj <^ 1, by |8C] 

(/!)3 



2M, 

-> oo, 

2M, 
oo, 



B, u;) 



(2Z + 1)!(2/)! 



(-4iMw) 



(4.157) 



(4.158) 



(4.159) 



If we insert these relations into the Green's functions, it is possible to carry out the integration in 
uj. In the near field limit 



1 



oj 



1 



167r2 

and in the far field limit 



(1 - 2M/r)tl 



, 1 V {I'-f /2M\ 
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oJ 
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WW vVJ 



OO. 



(4.160) 



(4.161) 



We immediately see that the Green's function is independent of the angular coordinates, as 
expected from the spherical symmetry. Note that the maximum value of / for which the expansion 
in Eqs. ( 4.16C| ) and ( 4.161 ) can be used depends upon the order of the leading terms which have 
been dropped in Eq. ( |4.159| ). If this correction is O ^(Mu;)''''^^ , then only the / = terms are 
significant, as Bq would then contain subdominant pieces which yield a contribution to GBEC^to) 
larger than the leading contribution from Bi. In what follows, we will explicitly retain only the 
I = contribution. In order to find the quantum inequality around a black hole we must evaluate 



Ap>--De GE{2to). 



(4.162) 



However, the only parts of the Euclidean box operator that are relevant are the temporal and radial 
terms, i.e., 

De^{1- 2M/r)-^dl + r-^dr[r^{l - 2M/r)dr]. (4.163) 

Upon taking the appropriate derivatives, and using the relation of the proper time of a stationary 
observer to the coordinate time, 

TO = toJl- 2M/r , (4.164) 



CHAPTER 4. SCALAR FIELD EXAMPLES 



65 



we find that the quantum inequahty is given by 



3 fl/2M\VTo\V 2M\-^ ( 2M\ \( 2M\ 



(4.165) 



as r — > 2M and 
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(4.166) 



as r — > oo. An alternative approach to finding the quantum inequality is to use the short time 



expansion from Section which yields 



A/5 > 



327r2r4 1671^2 



M2 



r4(l - 2M/r) 



+ 0(r2) + 



(4.167) 



Note that this short time expansion coincides with the first two terms of the r 2M form, 
Eq. (4.165). This is somewhat unexpected, as Eq. ( 4.165| ) is an expansion for small r — 2M with 
tq fixed, whereas Eq. ( [4.167| ) is an expansion for small tq with r fixed. 

We immediately see from Eq. ( 4.1661 ) that we recover the Minkowski space quantum inequality 
in the r ^ oo limit. If we consider experiments performed on the surface of the earth, where the 
radius of the earth is several orders of magnitude larger than its equivalent Schwarzschild radius, 
then the flat space inequality is an exceptionally good approximation. From Eq. ( |3.54| ), we can also 
find the proper sampling time for which the inequality Eq. ( [4.167] ) holds, 



To < 



2M 



r / 



(4.168) 



As was the case in two dimensions, if we allow the sampling time to go to infinity in the exact 
quantum inequality, we recover the QAWEC, Eq. ( 4.1501) , for the four-dimensional black hole. The 
QAWEC says that the renormalized energy density for an arbitrary particle state, sampled over 
the entirety of the rest observer's worldline, can never be more negative than the Boulware vacuum 
energy density. 



Chapter 5 

A Dynamic Spacetime: The Warp 
Drive 



In the preceding chapters, we looked at various static spacetimes in which quantum inequahties 
could be developed. Now we will apply the quantum inequality restrictions for a scalar field to 
Alcubierre's warp drive metric on a scale in which a local region of spacetime can be considered 
"flat." This is accomplished by taking the sampling time of the quantum inequality to be very 
small. Over short timescales, the equivalence principle tells us that the spacetime should be nearly 
static and almost flat. The short sampling time expansion of the quantum inequality developed in 
Section tells us that on short enough time scales we can apply the Minkowski space quantum 
inequality with reasonable reliability. Prom this we are able to place limits on the parameters of 
the "warp bubble." We will show that the bubble wall thickness is on the order of only a few 
hundred Planck lengths. Additionally, we will show that the total integrated energy density needed 
to maintain the warp metric with such thin walls is physically unattainable. 



5.1 Introduction 

In both the scientific community and pop culture, humans have been fascinated with the prospect 
of being able to travel between the stars within their own lifetimes. Within the framework of special 
relativity, the space-going traveler may move with any velocity up to, but not including, the speed 
of light. The astronauts would experience a time dilation which would allow them to make the 
round trip from earth to any star, and back to earth, in an arbitrarily short elapsed time. However, 
upon returning to earth such observers would find that their family and friends would have aged 



considerably. This is the well-known twin paradox P, 81 



Recently, Miguel Alcubierre proposed a metric [39|, fondly called the warp drive, in which a 
spaceship could travel to a star a distance D away and return home, such that the elapsed time for 
the stationary observers on earth would be less than 2D/c, where c is the velocity of light. What 
is most surprising about this spacetime is that the proper time of the space going traveler's trip is 
identical to the elapsed time on earth. However, the spaceship never locally travels faster than the 
speed of light. In fact, the spaceship can sit at rest with respect to the interior of the warp bubble. 
The ship is carried along by the spacetime, much in the same way that the galaxies are receding 
away from each other at extreme speeds due to the expansion of the universe, while locally they 
are at rest. The warp drive utilizes this type of expansion and contraction in order to achieve the 
ability to travel faster than light. 

Although warp drive sounds appealing, it does have one serious drawback. As with traversable 
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wormholes, in order to achieve warp drive we must employ exotic matter, that is, negative energy 
densities, which are a violation of the classical energy conditions. We have seen that the quantum 



inequality restrictions do allow negative energy densities to exist [16, 45, 47, 45]. However 



they place serious limitations on their magnitude and duration. Even before the short sampling 
time expansion of the quantum inequality was developed, researchers had applied the flat space 
quantum inequality to the curved spacetime geometries of wormholes |44] with the restriction that 
the negative energy be sampled on timescales smaller than the minimum local radius of curvature. 
It was argued that over such small sampling times, the spacetime would be locally flat and the 
inequalities would be valid. This led to the conclusion that static wormholes must either be on the 
order of several Planck lengths in size, or there would be large discrepancies in the length scales 
that characterize the wormhole. 

As we have seen in the preceding chapter, exact quantum inequalities can be developed for 
the static spacetimes in two, three, and four dimensions [^]. It was found that the quantum 
inequalities take the flat space form modified by a scale function which depends on the sampling 
time, the local radius of curvature, and the mass of the scalar field. In the limit of the sampling 
time being smaller than the local radius of curvature, the quantum inequalities reduce to the flat 
space form, often accompanied by higher order corrections due to the curvature |]35| , [3^ . In the 
limit of the radius of curvature going to infinity, the flat space inequalities are recovered. 

We would like to apply the same method to the warp drive metric, but such an exercise would 
require that we know the solutions to the Klein-Gordon equation for the mode functions of the scalar 
field. Such an approach, although exact, would be exceptionally difficult. In this section we will 
therefore apply the flat space inequality directly to the warp drive metric but restrict the sampling 
time to be small. We will thereby be able to show that the walls of the warp bubble must be 
exceedingly thin compared to its radius. This constrains the negative energy to an exceedingly thin 
band surrounding the spaceship. Similar results have been demonstrated for wormholes |^4|, where 
the negative energy is concentrated in a thin band around the throat. Recently, it has been shown 
for the Krasnikov metric |41], which also allows super luminal travel, that the required negative 
energy is also constrained to a very thin wall |42]. We will then calculate the total negative energy 
that would be required to generate a macroscopic sized bubble capable of transporting humans. As 
we will see, such a bubble would require physically unattainable energies. 



5.2 Warp Drive Basics 

Let us discuss some of the basic principles of the warp drive spacetime. We begin with a flat 
(Minkowski) spacetime and then consider a small spherical region, which we will call the bubble, 
inside this spacetime. On the forward edge of the bubble, we cause spacetime to contract, and on 
the trailing edge is an equal spacetime expansion. The region inside the bubble, which can be flat, 
is therefore transported forward with respect to distant objects. Objects at rest inside the bubble 
are transported forward with the bubble, even though they have no (or nominal) local velocity. 
Such a spacetime is described by the Alcubierre warp drive metric 

ds'^ = -dt^ + [dx - Vs{t)f{rs{t))dtf + V + (5-1) 

where Xs{t) is the trajectory of the center of the bubble and Vx{t) = dxs{t)/dt is the bubble's 
velocity. The variable rs{t) measures the distance outward from the center of the bubble given by 



rs{t) = ^(x-x,(t))2+y2+^2. 



(5.2) 
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The shape function of the bubble is given by f{rs), which Alcubierre originally chose to be 

tanh[a{rs - R)] - tanh[a{rs + R)] 
= 2tanh[aR] ' ^^'^^ 

The variable R is the radius of the warp bubble, and o" is a free parameter which can be used 
to describe the thickness of the bubble walls. In the large a limit, the function /(r^) quickly 
approaches that of a top hat function, where /(rg) = 1 for rg < R and zero everywhere else. It is 
not necessary to choose a particular form of /(r^). Any function will suffice so long as it has the 
value of approximately 1 inside some region < R and goes to zero rapidly outside the bubble, 
such that as Ts — > oo we recover Minkowski space. In order to make later calculations easier, we 
will use the piecewise continuous function 

1 r,<R-^, 
fpArs) = { -\{rs-R-^) R-^ <rs<R+^, (5.4) 

rs>R + ^, 

where R is the radius of the bubble. The variable A is the bubble wall thickness. It is chosen to 
relate to the parameter a for the Alcubierre form of the shape function by setting the slopes of the 
functions /(r^) and /p.c.('^s) to be equal at = R. This leads to 



l + tanh.^{aR) 



2 a ta.nh{aR) 



(5.5) 



which in the limit of large aR can be approximated by A ~ 2/ a. 

We now turn our attention to the solutions of the geodesic equation. It is straightforward to 
show that 

_ = n^ = (l,t,,(t)/(r,(t)),0,0), = (-1,0,0,0), (5.6) 

is a first integral of the geodesic equations. Observers with this four- velocity are called the Eulerian 
observers by Alcubierre. We see that the proper time and the coordinate time are the same for all 
observers. Also, the y and z components of the four-velocity are zero. The bubble therefore exerts 



no "force" in the directions perpendicular to the direction of travel. In Figure |5.1| , we have plotted 
one such trajectory for a observer that passes through the wall of a warp bubble at a distance p 
away from the center of the bubble. The x-component of the four- velocity is dependent on the shape 
function, and solving this explicitly for all cases can be rather difficult due to the time dependence 



of rs{t). A spacetime plot of an observer with the four- velocity given above is shown in Figure 5.2 
for a bubble with constant velocity. 

We see that the Eulerian observers are initially at rest. As the front wall of the bubble ap- 
proaches, the observer begins to accelerate in the direction of travel of the bubble, relative to 
observers at large distances. Once inside the bubble the observer moves with a nearly constant 
velocity given by 

dx{t) 



dt 



vs{tp)f{p), (5.7) 



which will always be less than the bubble's velocity unless p = (y^ -|- z^)^/^ = 0. The time tp is 
defined by rsitp) = p, i.e., it is the time at which the observer reaches the bubble equator. Such 
observers then decelerate, and are left at rest as they pass from the rear edge of the bubble wall. In 
other words no residual momentum is imparted to these observers during the "collision." However, 
they have been displaced forward in space along the trajectory of the bubble. 
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Figure 5.1: The path of an observer who passes through the outer region of the bubble, shown in 
the bubble's rest frame. As viewed from the interior of the bubble, the observer is moving to the 
left. 




Figure 5.2: The worldline (the dark line) of the geodesic observer passing through the outer region 
of a warp bubble, plotted in the observer's initial rest frame. The two lighter diagonal lines are the 
worldlines of the center of the bubble wall on the front and rear edges of the bubble, respectively. 
The bubble has a radius of 3, a velocity of 1, and a cr of 1. The plot shows an observer who begins 
at rest at x = 10, y^ + z^ = = 4. The shape function is of the form given by Alcubierre, Eq. ( ^.31 ). 
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There is also another interesting feature of these geodesies. As aheady noted, the observers will 
move with a nearly constant velocity through the interior of the bubble. This holds true for any 
value of p. However, the velocity is still dependent upon the value of p, so observers at different 
distances from the center of the bubble will be moving with different velocities relative to one 
another. If a spaceship of finite size is placed inside the bubble with its center of mass coincident 
with the center of the bubble, then the ship will experience a net "force" pushing it opposite to the 
direction of motion of the bubble so long as 1 — /(r<j) is nonzero at the walls of the ship. The ship 
would therefore have to use its engines to maintain its position inside the bubble. In addition, the 
ship would be subject to internal stresses on any parts that extended sufficiently far away from the 
center of the bubble 

In the above discussion we have used the Alcubierre form of the shape function, /(r^). If the 
piecewise continuous form, Eq. (^^), is used, one finds similar results with some modification. 
Inside the bubble, where < (R — A/2), every observer would move at exactly the speed of the 
bubble. So any observer who reached the bubble interior would continue on with it forever. This 
arises because everywhere inside the bubble, spacetime is perfectly flat due to /(r^) = 1. For 
observers whose geodesies pass solely through the bubble walls, so (i? — A/2) < p < {R + A/2), the 
result is more or less identical to that of the geodesies found with the Alcubierre shape function. 
This is the most interesting region because it contains the largest magnitude of negative energy. 

We now turn our attention to the energy density distribution of the warp drive metric. Using 
the first integral of the geodesic equations, it is easily shown that 

(T'^'^u u ) - (T'') - — G« - 1 ^s{t)p^ (dfirs)^ , . 

{I u,u,) -{I ) - - [ drs J ' ^ ^ ^ 

where p = [y^ + the radial distance perpendicular to the x-axis defined above. We 

immediately see that the energy density measured by any Eulerian observer is always negative, as 
shown in Alcubierre's original paper |3^. In Figure |5^ , we see that the distribution of negative 
energy is concentrated in a toroidal region perpendicular to the direction of travel. 



In Section 5A we will integrate the energy density over all of space to obtain the total negative 
energy required to maintain the bubble, under the restrictions of the quantum inequalities. As we 
will show, the total energy is physically unrealizable in the most extreme sense. 



5.3 Quantum Inequality Restrictions 

We begin with the quantum inequality (QI) for a free, massless scalar field in four-dimensional 



Minkowski spacetime originally derived by Ford and Roman [32| 

TO r {T^vu^'u^) 



vr J- 



oo 



^2 + ro2 ^"-"3^' ^'-'^ 
where r is an inertial observer's proper time, and tq is an arbitrary sampling time. It has been 



argued by Ford and Roman [44| that one may apply the QI to non-Minkowski spacetimes if the 
sampling time is smaller than the local radius of curvature. We know this to be true. It comes 
from the small sampling time expansion of the QI developed in Section 

We begin by taking the expression for the energy density ( |5.8D , and inserting it into the quantum 
inequality, Eq. ( ^.91) , yielding 

Vs{t)^ fdf{rs)\^ dt 



to / ^ 72—2 < ^ . (5.10) 

' 3 rj \ drs J V + to P % 
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Figure 5.3: The negative energy density is plotted for a longitudinal cross section of the warp metric 
traveling at constant velocity Vg = 1 to the right for the Alcubierre shape function. Black regions 
are devoid of matter, while white regions are maximal negative energy. 



If the time scale of the sampling is sufficiently small compared to the time scale over which the 
bubble's velocity is changing, then the warp bubble's velocity can be considered roughly constant, 
Vs{t) ~ Vb, during the sampling interval. We can now find the form of the geodesic at the time the 
sampling is taking place. Because of the small sampling time, the [t^ + ^o]~^ term becomes strongly 
peaked, causing the QI integral to sample only a small portion of the geodesic. We therefore place 
the observer at the equator of the warp bubble at t = 0. Then the geodesic is well-approximated 
by 

x{t) ~ fipHt , (5.11) 

which results in 

rsit)=[{vttf{f{p)-lf + pf^' . (5.12) 

Finally, we must specify the form of the shape function of the bubble. We can expand any shape 
function as a Taylor series about the sampling point, rs{t) — > p. Then we take the appropriate 
derivatives to obtain the term needed for the quantum inequality. We find 

dfirs 



dr. 



f'{p)+f"ip)[rs{t)-p] + ---. (5.13) 



The leading term is the slope of the shape function at the sampling point, which is in general 
roughly proportional to the inverse of the bubble wall thickness. We can therefore use, with no 



loss of generality, the piecewise continuous form of the shape function (5.4) to obtain a good order 
of magnitude approximation for any choice of shape function. The quantum inequality ( ^.10 ) then 
becomes 

'■+°° dt 3A2 



where 



Vh{l - fip)) 
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Formally the integral should not be taken over all time but just the time the observer is inside the 
bubble walls. However, the sampling function rapidly approaches zero. Therefore contributions 
to the integral from the distant past or the far future are negligible. The integral itself can be 
performed as the principal value of a contour that is closed in the upper half of the complex plane, 



/+0O 
-oo 



dt 



yielding an inequality of 



Vbto 



(1 - fip)) + 1 



(5.16) 



(5.17) 



The above inequality is valid only for sampling times on which the spacetime may be considered 
approximately flat. We must therefore find some characteristic length scale below which this occurs. 
For an observer passing through the bubble wall at a distance p from the center, the Riemann 
tensor in the static background frame can be calculated. Then we transform the components to 
the observer's frame by use of an orthonormal tetrad of unit vectors. In this frame, the tetrad is 
given by the velocity vector n''(t) and three unit vectors x, y, and z. The largest component of the 
Riemann tensor in the orthonormal frame is given by 



\R- 



This yields 



tyty I 



Trri.im. — 



1 



dp 
2A 



R 



tyty' 



V3 



Vb 



(5.18) 



(5.19) 



when y = p and the piecewise continuous form of the shape function is used. The sampling time 
must be smaller than this length scale, so we take 



to = a 



2A 



V3 



< a < 1. 



Vb 



(5.20) 



Here a is the ratio of the sampling time to the minimal radius of curvature. If we insert this into 
the quantum inequality and use 

^ ^ ^ « 1 , (5.21) 
P P 



we may neglect the term involving 1 — /(p) to find 



A < 



3 [3 Vb 



v2 ■ 



4 V TT 

Now as an example, if we let a = 1/10, then 

A < 10^ Vb Lpianck , 



(5.22) 



(5.23) 



where Lpianck is the Planck length. Thus, unless Vb is extremely large, the wall thickness cannot 
be much above the Planck scale. Typically, the walls of the warp bubble are so thin that the shape 
function could be considered a "step function" for most purposes. 
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5.4 Total Energy Calculation 

We will now look at the total amount of negative energy that is involved in the maintenance of 
a warp metric. For simplicity, let us take a bubble that moves with constant velocity such that 
Xs{t) = Vf)t. Because the total energy is constant, we can calculate it at time t = 0. We then have 

r,(t = 0) = [x^ + y^ + z^]^ =r. (5.24) 

With this in mind we can write the integral of the local matter energy density over proper volume 
as 

where g = Det\gij\ is the determinant of the spatial metric on the constant time hypersurfaces. 
Portions of this integration can be carried out by making a transformation to spherical coordinates. 
By doing so, we find 

1 2 r ifiiK^r' 



Since we are making only order of magnitude estimates of the total energy, we will use a piecewise 
continuous approximation to the shape function given by Eq. (5.4). Taking the derivative of this 



shape function, we find that the contributions to the energy come only from the bubble wall region, 
and we end up evaluating 

E = 




(5.27) 

For a macroscopically useful warp drive, we want the radius of the bubble to be at least in the 
range of 100 meters so that we may fit a ship inside. We showed in the previous section that the 
wall thickness is constrained by ( ^.23| ). If we use this constraint and let the bubble radius be equal 
to 100 meters, then we may neglect the second term on the right-hand-side of Eq. ( 5.27| ). It follows 
that 

E < -6.2 X 10™ Vb Lpianck ~ -6.2 X 10^^ Vb grams. (5.28) 
Because a typical galaxy has a mass of approximately 

MMiikyWay ~ 10^^ M^un = 2 X lO^^grams, (5.29) 
the energy required for a warp bubble is on the order of 

^ < -3 X 10^0 Mgalaxy Vb ■ (5.30) 

This is a fantastic amount of negative energy, roughly ten orders of magnitude greater than the 
total mass of the entire visible universe. 

If it were possible to violate the quantum inequality restrictions and make a bubble with a wall 
thickness on the order of a meter, things would be somewhat improved. The total energy required 
in the case of the same size radius and A = 1 meter would be on the order of a quarter of a solar 
mass, which would be more practical, yet still not attainable. 
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5.5 Summary 

We see from Eq. ( ^.231) , that quantum inequality restrictions on the warp drive metric constrain 
the bubble walls to be exceptionally thin. Typically, the walls are on the order of only hundreds 
or thousands of Planck lengths. Similar constraints on the size of the negative energy region have 
been found in the case of traversable wormholes . 

Prom Eq. ( |5.23| ), we may think that by making the velocity of the bubble very large we can make 
the walls thicker. However, this causes another problem. Por every order of magnitude by which 
the velocity increases, the total negative energy required to generate the warp drive metric also 
increases by the same magnitude. It is evident that for macroscopically sized bubbles to be useful 
for human transportation, even at subluminal speeds, the required negative energy is physically 
unattainable. 

On the other hand, we may consider the opposite regime. Warp bubbles are still conceivable 
if they are very tiny, i.e., much less than the size of an atom. Here the difference in length scales 
is not as great. As a result, a smaller amount of negative energy is required to maintain the warp 
bubble. Por example, a bubble with a radius of one electron Compton wavelength would require a 
negative energy of order E ~ —AOOMsun- 

The above derivation assumes that we are using a quantized, massless scalar field to generate 
the required negative energy. Similar quantum inequalities have been proven for both massive 
scalar fields ||3^, |4^ and the electromagnetic field [^Sj. In the case of the massive scalar field, the 
quantum inequality becomes even more restrictive, thereby requiring the bubble walls to be even 
thinner. Por the quantized electromagnetic field, the wall thickness can be made larger by a factor 
of \/2, due to the two spin degrees of freedom of the photon. However this is not much of an 
improvement over the scalar field case. 



Appendix A 

Inequalities 



In this appendix, we prove the following inequality: Let Aij be a real, symmetric n x n matrix 
with non-negative eigenvalues. (For the purposes of this manuscript, we may take either n = 2, 
for three-dimensional spacetimes, or n = 3, for four-dimensional spacetimes.) Further let h\ be 
a complex n-vector, which is also a function of the mode label A. Then in an arbitrary quantum 
state I'ip), the inequality states that 



A,A' 



Re Aij f^A*^i' (4"A'> ± h^xK, {axax')] > E ^u^A*M 



In order to prove this relation, we first note that 



V 



(A.l) 



(A.2) 



a=l 



where the Vj''"^ are the eigenvectors of A^j, and the Ka > are the corresponding eigenvalues. Now 
define the hermitian vector operator 



Note that 



Furthermore, 



Q' = J2 {hTai + h\ax) 



a=l 



(A.3) 
(A.4) 



a=l 



Q'^AijQA = 2ReY Aij \h\*h{,{a\ax>) + h\h{, {axax^)] + ^ Aijh\*h{ , (A.5) 



A, A' 



from which Eq. ( A.l ) with the '-|-'-sign follows immediately. The form of Eq. ( A.l ) with the '—'-sign 
can be obtained by letting /i-^ ih-'^. 

As a special case, we may take Aij = 5ij and obtain 



Rc^ [h*x-hy{a{ax')±hx-hy{ 



axax' 



A, A' 



> 



(A.6) 



As a further special case, we may take the vector to have only one component, e.g., hx 
{hx, 0, 0), in which case we obtain 



Re^ \hx*hx'{a{ax') ± hxhy {axay) > -^El^^l^- 



(A.7) 



A,A' 
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This last inequality was originally proven in [47| for real hx, and a simplified proof using the method 
adopted here is given in [48|. 



Appendix B 

Gamma Functions 



Here we discuss the following set of relations for gamma functions of complex argument, 



In addition we will prove the last two, Eqs. (|B.4| ) and ( p.Sl) . 

The first three gamma function relations, Eqs. ( [B.l - B.3| ), can be found in any mathematics text 
which includes special functions, such as Gradshteyn and Ryzhik |78], Section 8.33. Also, we note 
an important property for gamma functions of complex argument that is often not stated, 

r(x + iy)T{x — iy) = |r(3; + , for x and y real, (B.6) 

namely that the complex conjugate of a gamma function is equal to the gamma function of the 
complex conjugate of the argument. 

The last two gamma function relations, Eqs. ( |B.4| ) and ( p.5| ) are found from Eq. (8.332.4) of 
Gradshteyn and Ryzhik, which we repeat here: 

r(l + X + iy)T{l -x + iy)T{l + x- iy)T{l - x - iy) = ^ , (B.7) 

cosh lyn — cos zxvr 

for both X and y real. Using the relation for complex conjugates of Gamma functions we can rewrite 
the above expression as 

|r(l + x + ^yf = ^ , X l-'^-' + y') . (B.8) 

\T{l-x + iy)\^ cosh 2y7r- cos 2a:7r ^ ' 

In order to obtain relation ( |B.4D , we take x = 1/2 and y = k/2. We then have 

,2 _ 1 ^ 2vr^ [(1/2)^ + {k/ 2f 

|r(l/2 + i/c/2)|^ l + coshvr/;; 

Inserting Eq. (|B.2| ) into the above expression, we arrive at relation ( |B.4| ). Analogously, if we take 
X = 1 in Eq. ( [B.8D we obtain relation 



|r(3/2 + .A:/2)|2 = ——1-—-^ X '\ \ L7^:r'' ■ (B.9) 
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